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Tabulation of the Function «® = > 


sin ne 
n? 


Introduction. The function ¢(x) defined by 
* dt 
(1) (2) = fF togit +4 


occurs quite frequently in problems in Quantum Electrodynamics, as for example 
problems on vacuum polarization [1], scattering of light by light [2], and similar 
fourth order processes [3]. 

For real values of x, the function ¢(x) has been tabulated by Mitchell [4]. 
His tabulated function is ¢(0) — ¢(x), rather than our ¢(x) defined by (1). As 
in most literature he defines it as minus the integral from 0 to x. We shail however 
stick to the definition (1), which has already appeared in some of the literature [1 ]. 

The integrals appearing in the above mentioned problems cannot always be 
| expressed in terms of ¢(x) with the argument x real. In many cases there appear 
¢(x) with a complex x of absolute value 1. This happens when the integrals in 


1-6 i-3 
: ination of (75) and o(— 7-3) where 
question are expressed as combinations of ¢ its and ¢@ 14s where 6 


is to take all values from — © to + ©. When # is negative, the complex argument 
1 — 6/1 + 6 will have absolute value 1, and the problem reduces to that of find- 
ing ¢(e**). On use of the expansion 


-_ x (—)*t+1x* ¥ m 
@) axeet Sa Re 


the latter function can be expressed as 


(3) o(e*) = — #/4 + WO) 


where 


(4) v0) => ne 


1 


It was however shown by Kummer [5 ] that for complex x the imaginary part 
of ¢(x) can be expressed in terms of ¥(0) defined by equation (4). 

The function ¥(@), the Spence function, was tabulated by Clausen [6] for 
6 = 1, 2, ---, 180°. We find that these given intervals are not enough for easy 
computation of the function for any argument. 
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58 TABULATION OF A FUNCTION 


A tabuiation of the function ¥(6) for enough intervals will be useful for physi- 
cists to carry out explicit numerical evaluations of their calculations for compari- 
son with recent accurate experiment, which is necessary for the development of 
the theory. Therefore we found it worth while tabulating the function, given by 
Clausen for 180 intervals, for more intervals in such a way that the function can 
be easily computed for any value of the argument @. 


Some properties of the function (6). The series > ms 
1 





“ is slowly convergent 
n 


and the function is best computed by a numerical evaluation of the integral 


(5) v(@) = - in (2 cin Jat. 


The function is zero at @ = 0, @ = and has a maximum between these two 
limits at @ = 2/3. For @ < 2/2 the function can be expressed by the expansion 





(6) ¥@) = — Olog|o| +6+ 5 <—— 3, —" 
ne ~~ (n)! ~™" In(2n + 1)" 


For x/2 < @ < = the function can be expanded by use of the recurrence relation 


(7) v(x — 0) = ¥(0) — 49(28) 

and is given by 

8 (x — 6) = 61 2-y~ I"; ata th emanates 
(8) Vin = Olog 2 — 2. ony Ben" — Ye + 1) 


which shows that the function is linear in 6 for @ near x. 

Method of calculation. We start by evaluating the function from @ = 7 to 
6 = 2/2. The tabulated function logio (sin $¢) to 7 figures is used to obtain the 
value of the integrand in (5), at 540 intervals: 0, @:, ---. We have then used the 
known expression 


(9) Yn) = wL8F (00) +4061) + = + Ona) + 440.) — SF's) — FO] 


to compute ¥(6,). Here f(6;) is the integrand and w is the interval 1/1080. The 
remainder containing w* and higher powers was found not to affect the result 
expressed by 7 figures. 

The last tabulated value at @ = 2/2 was found to agree with the known 


) ; , 
; —— heck fi t of 
Catalan’s number, — (Qn + 1)” which was a good check on this first part o 


the work. 
The value of ¥(@) from @ = x/2 to @ = 0 will then be computed from the 
already calculated values on use of the relation ¥(6) = ¥(x — 0) + 4y(2@). On 
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reaching @ = 2/4 we had to start using newly calculated values of ¥(@) to complete 
the table until 6 = 0. For example the value at @ = 44° is ¥ (136°) + 4y(88°). 
The value at @ = 136° is contained in the first part of the work while the value at 
6 = 88° must be obtained from ¥(92°) + 4¥(176°), which is contained in the 
beginning of the second part of the work. 

The function was found to have the required maximum at @ = 2/3 and was 
further checked by directly evaluating it, using the expansion (6) for several 
values of @. 

The whole work has been carried out keeping the number of decimals 7, but 
it was found safer to tabulate it finally to 6 places of decimals which is enough 
for most purposes. 
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x 


¥ (x) 


0° 00’ 0.000000 


0.019897 
0.035761 
0.050103 
0.063457 
0.076075 
0.088108 
0.099654 
0.110783 
0.121547 
0.131988 
0.142137 
0.152022 
0.161663 
0.171081 
0.180291 
0.189307 
0.198141 
0.206803 
0.215304 
0.223652 
0.231855 
0.239919 
0.247851 
0.255656 
0.263340 
0.270908 
0.278364 
0.285712 
0.292956 
0.300100 
0.307147 
0.314101 
0.320963 
0.327738 
0.334427 
0.341032 
0.347558 
0.354004 
0.360374 
0.366670 
0.372893 
0.379045 
0.385128 
0.391143 
0.397092 
0.402977 
0.408799 
0.414558 
0.420258 
0.425898 
0.431480 
0.437005 
0.442474 
0.447888 
0.453249 
0.458557 
0.463813 
0.469018 
0.474173 
0.479279 


6 


19897 
15864 
14342 
13354 
12618 
12033 
11546 
11129 
10764 
10441 
10149 

9885 


9418 
9210 
9016 
8834 
8662 
850i 
8348 
8203 
8064 
7932 
7805 


7568 
7456 
7348 
7244 
7144 
7047 
6954 
6862 
6775 
6689 
6605 
6526 
6446 
6370 
6296 
6223 
6152 
6083 
6015 
5949 
5885 
5822 
5759 
5700 
5640 
5582 
5525 
5469 
5414 
5361 
5308 
5256 
5205 
5155 
5106 





x 


v(x) 


10° 00’ 0.479279 


10 
20 
30 
40 
50 
11°00 
10 
20 
30 
40 
50 
12°00 
10 
20 
30 
40 
50 
13°00 
10 
20 
30 
40 
50 
14°00 
10 
20 
30 
40 
50 
15°00 
10 
20 
30 
40 
50 
16° 00 
10 
20 
30 
40 
50 
17°00 
10 
20 
30 
40 
50 
18°00 
10 
20 
30 
40 
50 
19° 00 
10 
20 
30 
40 
50 
20° 00 


0.484336 
0.489346 
0.494309 
0.499227 
0.504098 
0.508925 
0.513709 
0.518449 
0.523146 
0.527801 
0.532415 
0.536988 
0.541521 
0.546014 
0.550468 
0.554883 
0.559261 
0.563600 
0.567903 
0.572168 
0.576398 
0.580592 
0.584750 
0.588874 
0.592963 
0.597018 
0.601040 
0.605028 
0.608983 
0.612906 
0.616797 
0.620656 
0.624483 
0.628280 
0.632045 
0.635781 
0.639486 
0.643161 
0.646807 
0.650424 
0.654012 
0.657571 
0.661102 
0.664605 
0.668080 
0.671527 
0.674948 
0.678341 
0.681708 
0.685048 
0.688362 
0.691650 
0.694912 
0.698149 
0.701361 
0.704547 
0.707709 
0.710846 
0.713958 
0.717047 
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6 


5057 
5010 
4963 
4918 
4871 
4827 
4784 
4740 
4697 
4655 
4614 
4573 
4533 
4493 
4454 
4415 
4378 
4339 
4303 
4265 
4230 
4194 
4158 
4124 
4089 
4055 
4022 
3988 
3955 
3923 
3891 
3859 
3827 
3797 
3765 
3736 
3705 
3675 
3646 
3617 
3588 
3559 
3531 
3503 
3475 
3447 
3421 
3393 
3367 
3340 
3314 
3288 
3262 
3237 
3212 
3186 
3162 
3137 
3112 
3089 





x 


10 
20 
30 
40 
50 
21°00 
10 
20 
30 
40 
50 
22° 00 
10 
20 
30 
40 
50 
23° 00 
10 
20 
30 
40 
50 
24° 00 
10 
20 
30 
40 
50 
25° 00 
10 
20 
30 
40 
50 
26° 00 
10 
20 
30 
40 
50 
27° 00 
10 
20 
30 
40 
50 
28° 00 
10 
20 
30 
40 
50 
29° 00 
10 
20 
30 
40 
50 
30° 00 





¥ (x) 
20° 00’ 0.717047 


0.720111 
0.723152 
0.726169 
0.729163 
0.732133 
0.735080 
0.738005 
0.740907 
0.743786 
0.746644 
0.749479 
0.752292 
0.755083 
0.757853 
0.760601 
0.763328 
0.766034 
0.768719 
0.771384 
0.774027 
0.776650 
0.779253 
0.781835 
0.784398 
0.786941 
0.789463 
0.791967 
0.794450 
0.796915 
0.799360 
0.801786 
0.804193 
0.806581 
0.808951 
0.811302 
0.813635 
0.815949 
0.818245 
0.820523 
0.822783 
0.825025 
0.827249 
0.829456 
0.831646 
0.833817 
0.835972 
0.838109 
0.840230 
0.842333 
0.844420 
0.846489 
0.848542 
0.850579 
0.852599 
0.854603 
0.856590 
0.858561 
0.860517 
0.862456 
0.864379 


6 


3041 
3017 
2994 
2970 
2947 
2925 
2902 
2879 
2858 
2835 
2813 
2791 
2770 
2748 
2727 
2706 
2685 
2665 
2643 
2623 
2603 
2582 
2563 
2543 
2522 
2504 
2483 
2465 
2445 
2426 
2407 
2388 
2370 
2351 
2333 
2314 
2296 
2278 
2260 
2242 
2224 
2207 
2190 
2171 
2155 
2137 
2121 
2103 
2087 
2069 
2053 
2037 
2020 
2004 
1987 
1971 
1956 
1939 
1923 








41 
17 


70 
47 
25 


79 
58 
35 
13 
91 


48 
27 


85 
65 
43 
123 
103 
582 
163 
34.3 
522 


183 
165 
145 
£26 
107 
388 
370 
351 
333 
314 
296 
278 
260 
242 
224 
207 


171 
155 
137 
121 
103 
087 
069 
053 
037 
020 
1004 
987 
971 
956 
939 
923 








x 


10 
20 


v(x) 
30° 00’ 0.864379 


0.866287 
0.868179 
0.870055 
0.871916 
0.873761 
0.875591 
0.877405 
0.879205 
0.880989 
0.882759 
0.884514 
0.886253 
0.887978 
0.889688 
0.891384 
0.893065 
0.894732 
0.896384 
0.898022 
0.899646 
0.901256 
0.902852 
0.904433 
0.906001 
0.907555 
0.909095 
0.910622 
0.912135 
0.913634 
0.915120 
0.916592 
0.918051 
0.919497 
0.920929 
0.922349 
0.923755 
0.925148 
0.926529 
0.927896 
0.929251 
0.930593 
0.931922 
0.933238 
0.934542 
0.935834 
0.937112 
0.938379 
0.939633 
0.940875 
0.942105 
0.943322 
0.944528 
0.945721 
0.946902 
0.948072 
0.949229 
0.950375 
0.951508 
0.952631 
0.953741 


1908 
1892 
1876 
1861 
1845 
1830 
1814 
1800 
1784 
1770 
1755 
1739 
1725 
1710 
1696 
1681 
1667 
1652 
1638 
1624 
1610 
1596 
1581 
1568 
1554 
1540 
1527 
1513 
1499 
1486 
1472 
1459 
1446 
1432 
1420 
1406 
1393 
1381 
1367 
1355 
1342 
1329 
1316 
1304 
1292 
1278 
1267 
1254 
1242 
1230 
1217 
1206 
1193 
1181 
1170 
1157 
1146 
1133 
1123 
1110 





x 


v(x) 


40° 00’ 0.953741 


10 
20 
30 
40 
50 
41°00 
10 
20 
30 
40 
50 
42°00 
10 
20 
30 
40 
50 
43° 00 
10 
20 
30 
40 
50 
44°00 
10 
20 
30 
40 
50 
45° 00 
10 
20 
30 
40 
50 
46° 00 
10 
20 
30 
40 
50 
47°00 
10 
20 
30 
40 
50 
48° 00 
10 
20 
30 
40 
50 
49° 00 
10 
20 
30 
40 
50 
50° 00 


0.954840 
0.955927 
0.957003 
0.958067 
0.959120 
0.960161 
0.961192 
0.962210 
0.963218 
0.964215 
0.965200 
0.966174 
0.967138 
0.968090 
0.969031 
0.969962 
0.970882 
0.971790 
0.972689 
0.973576 
0.974453 
0.975319 
0.976175 
0.977020 
0.977854 
0.978679 
0.979493 
0.980296 
0.981089 
0.981872 
0.982645 
0.983407 
0.984160 
0.984902 
0.985635 
0.986357 
0.987069 
0.987772 
0.988464 
0.989147 
0.989820 
0.990483 
0.991136 
0.991780 
0.992414 
0.993038 
0.993653 
0.994258 
0.994854 
0.995440 
0.996017 
0.996584 
0.997143 
0.997691 
0.998231 
0.998761 
0.999282 
0.999794 
1.000297 
1.000791 


61 


r) 


1099 
1087 
1076 


1053 
1041 
1031 
1018 
1008 
997 
985 
974 


952 
941 
931 
920 


899 
887 
877 





x 


v(x) 


50° 00’ 1.000791 


10 
20 
30 
40 
50 
51°00 
10 
20 
30 
40 
50 
52° 00 
10 
20 
30 
40 
50 
53°00 
10 
20 
30 
40 
50 
54°00 
10 
20 
30 
40 
50 
55° 00 
10 
20 
30 
40 
50 
56° 00 
10 
20 
30 
40 
50 
57° 00 
10 
20 
30 
40 
50 
58° 00 
10 
20 
30 
40 
50 
59° 00 
10 
20 
30 
40 
50 
60° 00 


1.001275 
1.001751 
1.002217 
1.002675 
1.003123 
1.003563 
1.003994 
1.004416 
1.004829 
1.005234 
1.005629 
1.006016 
1.006395 
1.006764 
1.007125 
1.007478 
1.007822 
1.008157 
1.008484 
1.008802 
1.009113 
1.009414 
1.009708 
1.009992 
1.010269 
1.010537 
1.010798 
1.011050 
1.011293 
1.011529 
1.011756 
1.011976 
1.012187 
1.012390 
1.012585 
1.012773 
1.012952 
1.013123 
1.013287 
1.013442 
1.013590 
1.013730 
1.013862 
1.013987 
1.014103 
1.014212 
1.014313 
1.014407 
1.014493 
1.014571 
1.014642 
1.014705 
1.014761 
1.014809 
1.014849 
1.014883 
1.014909 
1.014927 
1.014938 
1.014942 





v(x) 
1.014942 
1.014938 
1.014927 
1.014909 
1.014883 
1.014851 
1.014811 
1.014763 
1.014709 
1.014648 
1.014579 
1.014504 
1.014421 
1.014331 
1.014235 
1.014131 
1.014020 
1.013902 
1.013778 
1.013646 
1.013508 
1.013363 
1.013211 
1.013052 
1.012886 
1.012714 
1.012535 
1.012349 
1.012156 
1.011957 
1.011751 
1.011538 
1.011319 
1.011093 
1.010861 
1.010622 
1.010376 
1.010124 
1.009866 
1.009601 
1.009329 
1.009052 
1.008767 
1.008477 
1.008180 
1.007876 
1.007566 
1.007250 
1.006928 
1.006599 
1.006265 
1.005924 
1.005576 
1.005223 
1.004863 
1.004497 
1.004125 
1.003747 
1.003363 
1.002973 
1.002576 


8 


4 
11 
18 
26 
32 
40 
48 
54 
61 
69 
75 
83 
90 
96 

104 
111 
118 
124 
132 
138 
145 
152 
159 
166 
172 
179 
186 
193 


206 
213 
219 
226 
232 
239 
246 
252 
258 
265 





x v(x) 
70° 00’ 1.002576 
10 1.002174 
20 1.001766 
30 1.001351 
40 1.000931 
50 1.000505 
71°00 1.000072 
10 0.999634 
20 0.999190 
30 0.998740 
40 0.998284 
50 0.997822 
72°00 0.997355 
10 0.996881 
20 0.996402 
30 0.995917 
40 0.995426 
50 0.994930 
73°00 0.994427 
10 0.993919 
20 0.993406 
30 0.992886 
40 0.992361 
50 0.991831 
74°00 0.991294 
10 0.990752 
20 0.990205 
30 0.989652 
40 0.989093 
50 0.988529 
75°00 0.987959 
10 0.987384 
20 0.986803 
30 0.986217 
40 0.985625 
50 0.985028 
76°00 0.984425 
10 0.983817 
20 0.983204 
30 0.982585 
40 0.981961 
50 0.981331 
77°00 0.980696 
10 0.980056 
20 0.979411 
30 0.978760 
40 0.978104 
50 0.977443 
78°00 0.976776 
10 0.976104 
20 0.975427 
30 0.974745 
40 0.974058 
50 0.973366 
79°00 0.972668 
10 0.971965 
20 0.971257 
30 0.970544 
40 0.969826 
50 0.969103 
80°00 0.968375 
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402 
408 
415 
420 
426 
433 
438 


450 
456 
462 
467 
474 
479 
485 
491 
496 
503 
508 
513 
520 
525 
530 
537 
542 
547 
553 
559 
564 
570 
575 
581 
586 
592 
597 
603 
608 
613 
619 
624 
630 
635 
640 
645 
651 
656 
661 
667 
672 
677 
682 
687 
692 
698 
703 
708 
713 
718 
723 
728 








x v(x) 
80° 00’ 0.968375 
10 0.967642 
20 0.966903 
30 0.966160 
40 0.965412 
50 0.964658 
81°00 0.963900 
10 0.963137 
20 0.962369 
30 0.961596 
40 0.960818 
50 0.960035 
82°00 0.959247 
10 0.958455 
20 0.957657 
30 0.956855 
40 0.956048 
50 0.955236 
83°00 0.954419 
10 0.953598 
20 0.952771 
30 0.951940 
40 0.951105 
50 0.950264 
84°00 0.949419 
10 0.948569 
20 0.947714 
30 0.946855 
40 0.945991 
50 0.945123 
85°00 0.944249 
10 0.943372 
20 0.942489 
30 0.941602 
40 0.940711 
50 0.939814 
86°00 0.938914 
10 0.938008 


20 0.937099 - 


30 0.936184 
40 0.935266 
50 0.934342 
87°00 0.933415 
10 0.932482 
20 0.931546 
50 0.930605 
40 0.929659 
50 0.928709 
88°00 0.927755 
10 0.926796 
20 0.925833 
30 0.924866 
40 0.923894 
50 0.922918 
89°00 0.921938 
10 0.920953 
20 0.919964 
30 0.918971 
40 0.917973 
50 0.916972 
90°00 0.915966 


733 
739 
743 
748 
754 
758 
763 
768 
773 
778 
783 
788 
792 
798 
802 
807 
812 
817 
821 
827 
831 
835 
841 
845 
850 
855 
859 
864 
868 
874 
877 
883 
887 
891 
897 


915 


924 
927 


1001 
1006 








x 





v(x) 


90° 00 0.915966 
733 10 0.914955 
739 20 0.913941 
743 30 0.912922 
748 40 0.911899 
754 50 0.910872 
758 91°00 0.909841 
763 10 0.908806 
768 20 0.907766 
730 30 0.906722 
778 40 0.905675 
783 50 0.904623 
788 92°00 0.903567 
792 10 0.902507 
798 20 0.901443 
802 30 0.900374 
807 40 0.899302 
812 50 0.898226 
817 93°00 0.897145 
821 10 0.896061 
827 20 0.894973 
831 30 0.893880 
835 40 0.892784 
841 50 0.891684 
845 # 94°00 0.890580 
850 10 0.889471 
855 20 0.888359 
859 30 0.887243 
864 40 0.886123 
868 50 0.884999 
874 95°00 0.883872 
877 10 0.882740 
883 20 0.881605 
887 30 0.880465 
891 40 0.879322 
897 50 0.878175 
900 96°00 0.877024 
906 10 0.875870 
909 20 0.874711 
915 30 0.873549 
918 40 0.872383 
924 50 0.871213 
927 97°00 0.870040 
933 10 0.868862 
936 20 0.867681 
941 30 0.866497 
946 40 0.865308 
950 50 0.864116 
954 98°00 0.862920 
959 10 0.861721 
963 20 0.860518 
967 30 0.859311 
972 40 0.858100 
976 50 0.856886 
980 99°00 0.855668 
985 10 0.854447 
989 20 0.853222 
993 30 0.851994 
998 40 0.850762 
1001 50 0.849526 
1006 100° 00 0.848287 








1104 
1109 
1112 
1116 
1120 
1124 
1127 
1132 
1135 
1140 
1143 
1147 
1151 
1154 
1159 
1162 
1166 
1170 
1173 
1178 
1181 
1184 
1189 
1192 
1196 
1199 
1203 
1207 
1211 
1214 
1218 
1221 
1225 
1228 
1232 
1236 
1239 





x ¥ (x) 
100° 00’ 0.848287 
10 0.847044 
20 0.845797 
30 0.844548 
40 0.843294 
50 0.842037 
101°00 0.840777 
10 0.839513 
20 0.838246 
30 0.836975 
40 0.835701 
50 0.834423 
102°00 0.833142 
10 0.831857 
20 0.830569 
30 0.829278 
40 0.827983 
50 0.826685 
103° 00 0.825383 
10 0.824078 
20 0.822770 
30 0.821458 
40 0.820143 
50 0.818825 
104°00 0.817503 
10 0.816179 
20 0.814850 
30 9.813519 
40 0.812184 
50 0.810846 
105°00 0.809505 
10 0.808160 
20 0.806813 
30 0.805462 
40 0.804107 
50 0.802750 
106° 00 0.801389 
10 0.800025 
20 0.798658 
30 0.797288 
40 0.795915 
50 0.794538 
107° 00 0.793159 
10 0.791776 
20 0.790390 
30 0.789001 
40 0.787609 
50 0.786214 
108°00 0.784816 
10 0.783414 
20 0.782010 
30 0.780603 
40 0.779192 
50 0.777779 
109° 00 0.776362 
10 0.774943 
20 0.773520 
30 0.772094 
40 0.770666 
50 0.769234 
110°00 0.767800 
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1243 
1247 
1249 
1254 
1257 
1260 
1264 
1267 
1271 
1274 
1278 
1281 
1285 
1288 
1291 
1295 
1298 
1302 
1305 
1308 
1312 
1315 
1318 
1322 
1324 
1329 
1331 
1335 
1338 
1341 
1345 
1347 
1351 
1355 
1357 
1361 
1364 
1367 
1370 
1373 
1377 
1379 
1383 
1386 
1389 
1392 
1395 
1398 
1402 
1404 
1407 
1411 
1413 
1417 
1419 
1423 
1426 
1428 
1432 
1434 





x 


¥(x) 


110° 00’ 0.767800 


10 
20 


0.766362 
0.764922 
0.763478 
0.762032 
0.760583 
0.759131 
0.757676 
0.756218 
0.754757 
0.753293 
0.751826 
0.750357 
0.748885 
0.747410 
0.745932 
0.744451 
0.742967 
0.741481 
0.739991 
0.738499 
0.737004 
0.735507 
0.734006 
0.732503 
0.730997 
0.729489 
0.727978 
0.726463 
0.724947 
0.723427 
0.721905 
0.720380 
0.718853 
0.717322 
0.715790 
0.714254 
0.712716 
0.711175 
0.709632 
0.708086 
0.706537 
0.704986 
0.703432 
0.701875 
0.700316 
0.698755 
0.697191 
0.695624 
0.694055 
0.692483 
0.690909 
0.689332 
0.687753 
0.686171 
0.684586 
0.683000 
0.681410 
0.679819 
0.678224 
0.676628 





—$ 


1438 
1440 
1444 
1446 
1449 
1452 
1455 
1458 
1461 
1464 
1467 
1469 
1472 
1475 
1478 
1481 
1484 
1486 
1490 
1492 
1495 
1497 
1501 
1503 
1506 
1508 
1511 
1515 
1516 
1520 
1522 
1525 
1527 
1531 
1532 
1536 
1538 
1541 
1543 
1546 
1549 
1551 
1554 
1557 
1559 
1561 
1564 
1567 
1569 
1572 
1574 
1577 
1579 
1582 
1585 
1586 
1590 
1591 
1595 
1596 





x v(x) 
120° 00’ 0.676628 
10 0.675029 
20 0.673427 
30 0.671823 
40 0.670217 
50 0.668608 
121°00 0.666997 
10 0.665383 
20 0.663767 
30 0.662149 
40 0.660528 
50 0.658905 
122°00 0.657280 
10 0.655652 
20 0.654022 
30 0.652390 
40 0.650755 
50 0.649118 
123°00 0.647478 
10 0.645837 
20 0.644193 
30 0.642547 
40 0.640898 
50 0.639247 
124°00 0.637594 
10 0.635939 
20 0.634281 
30 0.632622 
40 0.630960 
50 0.629295 
125°00 0.627629 
10 0.625960 
20 0.624290 
30 0.622617 
40 0.620941 
50 0.619264 
126°00 0.617585 
10 0.615903 
20 0.614219 
30 0.612533 
40 0.610845 
50 0.609155 
127°00 0.607463 
10 0.605768 
20 0.604072 
30 0.602373 
40 0.600672 
50 0.598969 
128°00 0.597265 
10 0.595558 
20 0.593849 
30 0.592138 
40 0.590424 
50 0.588709 
129°00 0.586992 
10 0.585273 
20 0.583552 
30 0.581829 
40 0.580104 
50 0.578376 
130°00 0.576647 
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1599 
1602 
1604 
1606 
1609 
1611 
1614 
1616 
1618 
1621 
1623 
1625 
1628 
1630 
1632 
1635 
1637 
1640 
1641 
1644 
1646 
1649 
1651 
1653 
1655 
1658 
1659 
1662 
1665 
1666 
1669 
1670 
1673 
1676 
1677 
1679 
1682 
1684 
1686 
1688 
1690 
1692 
1695 
1696 
1699 
1701 
1703 
1704 
1707 
1709 
1711 
1714 
1715 
1717 
1719 
1721 
1723 
1725 
1728 
1729 





x ¥ (x) 
130° 00’ 0.576647 
10 0.574916 
20 0.573183 
30 0.571448 
40 0.569711 
50 0.567972 
131°00 0.566231 
10 0.564489 
20 0.562744 
30 0.560997 
40 0.559249 
50 0.557498 
132°00 0.555746 
10 0.553992 
20 0.552236 
30 0.550478 
40 0.548718 
50 0.546956 
133°00 0.545193 
10 0.543427 
20 0.541660 
30 0.539891 
40 0.538120 
50 0.536348 
134°00 0.534573 
10 0.532797 
20 0.531019 
30 0.529239 
40 0.527458 
50 0.525674 
135°00 0.523889 
10 0.522102 
20 0.520314 
30 0.518524 
40 0.516731 
50 0.514938 
136°00 0.513142 
10 0.511345 
20 0.509546 
30 0.507745 
40 0.505943 
50 0.504139 
137°00 0.502334 
10 0.500526 
20 0.498717 
30 0.496907 
40 0.495094 
50 0.493280 
138°00 0.491465 
10 0.489648 
20 0.487829 
30 0.486008 
40 0.484186 
50 0.482363 
139°00 0.480538 
10 0.478711 
20 0.476883 
30 0.475053 
40 0.473221 
50 0.471388 
140°00 0.469554 
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1731 
1733 
1735 
1737 
1739 
1741 
1742 
1745 
1747 
1748 
1751 
1752 
1754 
1756 
1758 
1760 
1762 
1763 
1766 
1767 
1769 
1771 
1772 
1775 
1776 
1778 
1780 
1781 
1784 
1785 
1787 
1788 
1790, 
1793 
1793 
1796 
1797 
1799 
1801 
1802 
1804 
1805 
1808 
1809 
1810 
1813 
1814 
1815 
1817 
1819 
1821 
1822 
1823 
1825 
1827 
1828 
1830 
1832 
1833 
1834 








x v(x) 
140° 00’ 0.469554 
10 0.467717 
20 0.465880 
30 0.464041 
40 0.462200 
50 0.460358 
141°00 0.458514 
10 0.456669 
20 0.454822 
30 0.452974 
40 0.451124 
50 0.449273 
142°00 0.447421 
10 0.445567 
20 0.443711 
30 0.441854 
40 0.439996 
50 0.438136 
143°00 0.436275 
10 0.434412 
20 0.432548 
30 0.430683 
40 0.428816 
50 0.426947 
144°00 0.425078 
10 0.423207 
20 0.421334 
30 0.419461 
40 0.417586 
50 0.415709 
145°00 0.413831 
10 0.411952 
20 0.410072 
30 0.408190 
40 0.406307 
50 0.404423 
146°00 0.402537 
10 0.400650 
20 0.398762 
30 0.396872 
40 0.394982 
50 0.393090 
147°00 0.391196 
10 0.389302 
20 0.387406 
30 0.385509 
40 0.383610 
50 0.381711 
148°00 0.379810 
10 0.377908 
20 0.376005 
30 0.374101 
40 0.372195 
50 0.370288 
149°00 0.368380 
10 0.366471 
20 0.364561 
30 0.362649 
40 0.360737 
50 0.358823 
150°00 0.356908 


—3 


1837 
1837 
1839 
1841 
1842 


1845 
1847 
1848 
1850 
1851 
1852 
1854 
1856 
1857 
1858 
1860 
1861 
1863 
1864 
1865 
1867 
1869 
1869 
1871 
1873 
1873 
1875 
1877 
1878 
1879 
1880 
1882 
1883 
1884 
1886 
1887 
1888 


“1890 


1890 
1892 
1894 
1894 
1896 
1897 
1899 
1899 
1901 
1902 
1903 
1904 
1906 
1907 
1908 
1909 
1910 
1912 
1912 
1914 
1915 





1! 


1! 


1: 


1 


1! 


i 











x 


10 
20 
30 
40 
50 
151° 00 
10 
20 


¥(x) 
150° 00’ 0.356908 


0.354992 
0.353075 
0.351157 
0.349238 
0.347317 
0.345396 
0.343473 
0.341549 
0.339624 
0.337698 
0.335771 
0.333843 
0.331914 
0.329984 
0.328053 
0.326121 
0.324187 
0.322253 
0.320318 
0.318382 
0.316444 
0.314506 
0.312567 
0.310626 
0.308685 
0.306743 
0.304800 
0.302856 
0.300911 
0.298965 
0.297018 
0.295070 
0.293121 
0.291171 
0.289220 
0.287269 
0.285316 
0.283363 
0.281409 
0.279454 


50 0.277498 


0.275541 
0.273583 
0.271625 
0.269665 
0.267705 
0.265744 
0.263782 
0.261819 
0.259856 
0.257891 
0.255926 
0.253960 
0.251993 
0.250026 
0.248057 
0.246088 
0.244118 
0.242148 
0.240176 
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1916 
1917 
1918 
1919 
1921 
1921 
1923 
1924 
1925 
1926 
1927 
1928 
1929 
1930 
1931 
1932 
1934 
1934 
1935 
1936 
1938 
1938 
1939 
1941 
1941 
1942 
1943 
1944 
1945 
1946 
1947 
1948 
1949 
1950 
1951 
1951 
1953 
1953 
1954 
1955 
1956 
1957 
1958 
1958 
1960 
1960 
1961 
1962 
1963 
1963 
1965 
1965 
1966 
1967 
1967 
1969 
1969 
1970 
1970 
1972 





x 


v(x) 


160° 00’ 0.240176 


10 
20 
30 
40 
50 
161° 00 
10 
20 
30 
40 
50 
162° 00 
10 
20 
30 
40 
50 
163° 00 
10 
20 
30 
40 
50 
164° 00 
10 
20 
30 
40 
50 
165° 00 
10 
20 
30 
40 
50 
166° 00 
10 
20 
30 
40 
50 
167° 00 
10 
20 
30 
40 
50 
168° 00 
10 
20 
30 
40 
50 
169° 00 
10 
20 
30 
40 
50 
170° 00 


0.238204 
0.236231 
0.234258 
0.232283 
0.230308 
0.228333 
0.226356 
0.224379 
0.222401 
0.220423 
0.218443 
0.216463 
0.214483 
0.212502 
0.210520 
0.208537 
0.206554 
0.204570 
0.202586 
0.200601 
0.198615 
0.196629 
0.194642 
0.192654 
0.190666 
0.188677 
0.186688 
0.184698 
0.182708 
0.180717 
0.178725 
0.176733 
0.174740 
0.172747 
0.170753 
0.168759 
0.166764 
0.164769 
0.162773 
0.160777 
0.158780 
0.156783 
0.154785 
0.152787 
0.150788 
0.148789 
0.146789 
0.144789 
0.142789 
0.140788 
0.138786 
0.136784 
0.134782 
0.132780 
0.130777 
0.128773 
0.126769 
0.124765 
0.122760 
0.120755 
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x v(x) 
170° 00’ 0.120755 
10 0.118750 
20 0.116744 
30 0.114738 
40 0.112732 
50 0.110725 
171°00 0.108718 
10 0.106710 
20 0.104702 
30 0.102694 
40 0.100686 
50 0.098677 
172°00 0.096668 
10 0.094659 
20 0.092649 
30 0.090639 
40 0.088629 
50 0.086619 
173° 00 0.084608 
10 0.082597 
20 0.080586 
30 0.078574 
40 0.076562 
50 0.074551 
174°00 0.072538 
10 0.070526 
20 0.068513 
30 0.066500 
40 0.064487 
50 0.062474 
175°00 0.060461 
10 0.058447 
20 0.056433 
30 0.054419 
40 0.052405 
50 0.050391 
176° 00 0.048377 
10 0.046362 
20 0.044347 
30 0.042332 
40 0.040317 
50 0.038302 
177°00 0.036287 
10 0.034272 
20 0.032256 
30 0.030241 
40 0.028225 
50 0.026209 
178°00 0.024194 
10 0.022178 
20 0.020162 
30 0.018146 
40 0.016130 
50 0.014114 
179°00 0.012097 
10 0.010081 
20 0.008065 
30 0.006049 
40 0.004033 
50 0.002016 
180° 00 0.000000 


—8 


2005 
2006 
2006 
2006 
2007 
2007 
2008 
2008 
2008 
2008 
2009 
2009 
2009 
2010 
2010 
2010 
2010 
2011 
2011 
2011 
2012 
2012 
2011 
2013 
2012 
2013 
2013 
2013 
2013 
2013 
2014 
2014 
2014 
2014 
2014 
2014 
2015 
2015 
2015 
2015 
2015 
2015 
2015 
2016 
2015 
2016 
2016 
2015 
2016 
2016 
2016 
2016 
2016 
2017 
2016 
2016 
2016 
2016 
2016 
2016 
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A Modified Chebyshev-Everett Interpolation 
Formula 


1. Introduction. It is a well-known property of the Chebyshev polynomials 
that, of all polynomials with leading coefficient unity, they possess the smallest 
absolute upper bound when the argument is allowed to vary between their limits 
of orthogonality. This property suggests the use of such polynomials as a means 
of interpolation. As is explained in Kopal’s recent book on numerical analysis [1], 
a power series rearranged as a series of Chebyshev polynomials is reduced to an 
economic form, insofar as the “‘maximum”’ accuracy is attained with a very small 
number of terms. We use Kopal’s notation [1] here. 

The purpose of this note is to rearrange the Everett interpolation formula, 


f(m) = f(O) + my + Eotidcti + EyitAyii + --- + Ey @mAg@m) 4 By @a)Q, @n) 


as a series of suitable Chebyshev polynomials, thereby reducing considerably the 
number of terms depending on m, required to give a prescribed accuracy, while 
retaining the use of even order differences only. Following Miller [2] the 
polynomials C2;,:(2 — 2m) and C2;4:(2m) could be used, where C;(2m) 
= 2 cos (j cos m). A disadvantage in doing so is the non-vanishing of such 
polynomials at tabular points. To overcome this difficulty Kopal [1] and Miller 
[7] have suggested an expansion in terms of the polynomials 


TT 


4j +2 








Coj41 | (2 — 2m)cos and Cassa {2m cos / . 


4j +2 

Such an expansion Kopal calls the Modified Chebyshev-Everett Interpolation 
Formula and states that “--- it represents potentially the most powerful and 
economic interpolation formula which can as yet be devised” [1]. See also 
Clenshaw and Olver [9]. 

: Explicitly, following [1], p. 513, let 





(1) f(m) = f(0) + mA; + z= [soiCans {2 — 2m) cos rr * | 


j=l 





+ B19Coj41 {2m cos rr : | 


secti+ ( . ) 
" aj +2 














| + 19 Caj41|2m cos ~ | 


x| 18a (2 — 2m)cos +2 


T 
4j +2 


where the coefficients No? and N,°* will be designated modified differences of 
the 2j-th order. Owing to the symmetry of the Everett coefficients E®*) and to the 
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fact that C2;4:(x) is an odd function of x, the N,°@*’s will be the same function of 
the A1’s as the No*’s are of the Ay’s. 

2. Calculation of the Modified Differences N,,,°°”. The coefficients No, ;°” 
are determined from the equivalence of the two expressions 


f(m) = f(0) + mA; + > (Eo®? Ag?) + E, 24,0) 


i=l 


2j+1 = ) 
a Gas 


= f(0) + md; + ~ PQ; + 1)! 














x| 19a (2 — 2m)cos 


ncnleeegtal 
ij sl 4 N,i PCays|2m cos 


27+1 bd ) 
—— (G5 


3 E,@),a,¢@) = —— 
( ) = 1 1 ~ 224+1 (27 oa 1)! 


or, in effect, from 











Ms°Caya| 2m cos 4j : | ° 
Both sides of (3) are expanded in powers of m and the coefficients of equal powers 
equated. The procedure is quite straightforward, but the algebra is rather un- 
wieldy, and is not reproduced here. 

The following general relations can be obtained by equating the coefficients 
of m**-* and m: 


1jn 1 T 
(Qa—4) — (2a—4) — — os dtp 2 (22—2) 
” . i(: 2(n —1) ai5)4 














is s 
+ + nGn +0 el nt? "in—2 











1440 ~ 64n(2n — 1) 32(m — 1) 
ma+i 1 o£ ) " std 
x( 3. mn in+2 } ao» 
and 
237 
<a Gy. ane ;___ 9B t2 vy an 
Sai pie x (Oa 


These relations are of assistance in computing the modified differences. 
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If differences up to the tenth are retained the following closed forms for the 
coefficients Nii, Niv, ---, N* are obtained: 


N= = 4+, 
Nviii = Aviii + a;A*, 























Nvi = Avi + a,Avili + 5A, 
Nv = Aiv a a;Avi a bAviii + c,A%, 
Nii = Aii oo a,Ai¥ -b b3Avi - cy Aviii ao d;A*, 
where 
a;>- 1 (i+? om 1. 2 bad ) 
' 4\ 3 3” G+2)' 
oe sect ball 
», - GtnGi+e “Fa, WG-2 
P 1440 64j(27-1) * 8(j-—1)' 
139 1 T 3 T bs T 
oe oe eee ccmmenenaiemnattiin at ens o ian aa ie 
*" 6048) 5-3-9 99 — 7.985 7g F 9g 8° Ty 
41 5 T a Tr bs T 
ae eS ee eee Ie ee ee 
os 3024 * 16-81% 18 960° 14 * 16 10” 
479 15 T 5das5 T bs T Cs T 
d= _ Sie 6 4 veg 
* = 751200 27-101°** 22> + 


16-8! 18 8-5!" 14. 16° 10 


The numerical values of the coefficients a, b, c, d-are given in Table I. 


2j7+1 = ) 
—- (=, 


22127 +1)!’ 
Chebyshev polynomials §o,:‘” in equation (1) are obtained. These are given in 
Table II (4 = 0,1). 

As will be seen, the coefficients 8‘) diminish very rapidly and in practice it 
will only be necessary to use about two terms of the series (1). Should it be found 
unnecessary to retain the higher differences, they should simply be omitted from 
the 8 or N@?, 

3. Besselian and Comrie-type Modified Differences. Clearly, the well-known 
Bessel interpolation formula can also be regrouped in terms of the Chebyshev 
polynomials, and the modified differences corresponding to it obtained in a similar 
manner. Modified differences arising in connection with Bessel’s formula were 
first introduced into computational practice by E. W. Brown [3], and subse- 
quently (but independently) by Camp [4] and Comrie [5]. Comrie, in particular, 
has done so much to propagate their use that they are usually associated with his 
name. It must be emphasized that the Comrie type of (Besselian) modified 
differences is not identical with what may be called the Chebyshev-Bessel modi- 
fied differences, which would be obtained by an analogous procedure to that 
developed in the preceding sections of this note. 
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Comrie used modified differences M‘ whose coefficients were determined 
to minimize the error of the Bessel interpolation formula truncated after a given 
number of terms. Comrie never published details of the process by which these 
constants were determined ; his argument has been reconstructed in [1], Chapter 
II. Although these constants were determined by Comrie with specific reference 
to Bessel’s formula, it has become customary to use them in Everett’s formula 
also. In this case an expansion 


f(m) = f(O) + mA; + EotiAot + E,4A;4 + --- 
+ Eo?) (Ap? + AjAo@#® + BjAp@i +4 ---) 
+ E,@)(A,@) 4+ A,;A,e#®) + BA, CH + «..) 


is used, where the coefficients A;, B; are those determined by Comrie for the 


Bessel formula. 
For convenience, write 


27+1 ba ) 
= (3 2 
2i(j + 1)! 











Coj41 (2m cos 


T2j+1 = 


+3) 
4j+2 

Consider a specific example where differences up to the eighth are retained. 
Then f(m) may be expanded: 
(4)  f(m) = f(0) + mA, + Ei%Ai4 + EyvAy + Evia! + Eyrittayviti 


+ similar terms in EoAo, 


(5) f(m) = f(0) + mA; + r3(A,ii + adiv + BAI + cAyvii) 
+ 75(Ari¥ + dAy! + eA;"ii) 
~ 7 (Ay¥i + fAyviii) 
+ rpAyvili + oe, 


A comparison of these formulae shows that 


rs = Ei 


ar; + 75 = Ei’. 
Now if we regard the last three terms of (5) as error terms then 


f(m) ~ f(0) + mA, + rs(Arti + aAriv + BA! + cAvitt) + --. 
= f(0) + mA; + Extn, + ---, 


Reference to Comrie’s paper [5] shows that N;" is in fact identical with the 
modified second difference used by Comrie. 

However, there seems to be no simple relation between N;' and higher order 
modified differences, and the higher Comrie modified differences. 
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4. Example. As an example, an interpolation has been made in a table of 
Fermi-Dirac functions at present being prepared by the author. 

A portion of this table, together with ist, 2nd, 4th, 6th, 8th, and 10th differ- 
ences, is reproduced below. It has been chosen because the differences do not fall 
off rapidly, as can be seen. 


n F(n) Ay Al Av avi Avil at 
—2.0 0.1145 8783 
6921 403 
—1.5 1838 0186 3748 501 297 363 —135 029 65 578 —20 496 
10669 904 
-1.0 2905 0090 5259 311 130 749  —125 114 110 942  -—119 432 
15929 215 
—0.5 4497 9305 6900 870 -—160 979 —4 256 36 873 —68 774 
22830 085 
0.0 6780 9389 asi) sae 8381 450 -—456 963 153 475 -—105 970 113 458 
0.5 9902 0925 9405 068  -—599 471 205 237 —135 355 132 585 
40616 603 
1.0 1.3963 7528 9829 215 —536 743 121 643 —32 156 —17 467 
50445 818 
1.5 1.9008 3346 9716 618 —352 371 5 895 53 577 —82 922 
60162 437 


2.0 2.5024 5783 
From this table we find, e.g., the value of F(—0.3) from the formula: (m=0.4), 
F(—0.3) = F(—0.5) + 0.4A, + Eo't(0.4) Nol + Ey#(0.4) Ni 


2 =) 2 ( =) 
+ Bo cx(1.2 cos 10 + BC; | 0.8 cos 10 
where 


Not = 0.06929 411 Bo? = — 0.00000 052 
N,it = 0.08470 882 B62 = — 0.00000 171. 


The values of Cs (12 cos =) and Cs (os cos =) are obtained by linear inter- 


polation in the table of Chebyshev polynomials [6]. 

For comparison, columns (1) and (2) below give respectively the value of 
F(—0.3), as calculated with the same number of terms in the interpolation 
formula, by the straight Everett formula, 


f(m) = f(0) oo may EgtiAgit a E,'iA, ii os Eo'* Aoi*¥ — E,'A,i* 


and the modified Everett formula using Comrie’s coefficients in the modified 
fourth difference, 


f(m) = f(0) + md, + Eo'tAc! + E,%Ai" + Eo Mo + Ey" My 


where Mi, = Aly, — 0.207A5;, + 0.045457. Column (3) gives the value com- 
puted by the present method. The value obtained by direct calculation is 
0.53193 157. It will be seen that the present method gives substantially better 
results for this case. 
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1. 2. 3. 
F(—0.3) =  0.44979305 0.44979305 0.44979305 
+ 9132034 9132034 9132034 
- 441656 441656 443482 
- 469361 469361 474369 
- 1875 1845 11 
- 4913 5306 318 





0.53193534 0.53193171 0.53193159 





In most applications high order differences will fall off rapidly and probably 
only the modified second difference will be required. If this is calculated and 
printed with a table, interpolation becomes a very simple procedure. 

5. Use of the Expansion. An interesting discussion about the merits of a 
Chebyshev-Everett interpolation formula is recorded in Appendix 11 of reference 
[7] to which the reader is referred. (This report was received after the major part 
of this note had been written.) 

It should be emphasized that the chief advantage of this type of expansion is 
that the number of terms in the interpolating expansion of a function f(m), 
dependent on m, is reduced considerably. Clearly, the expansion will have its 
greatest use when several interpolations within a given interval are required. 
The complicated modified differences, once calculated, are availabie for all inter- 
polations within the relevant interval, the number of terms dependent on m 
being usually only two, viz., 


Cs {2m cos = | and C3} (2 — 2m) cos 7 


A table of the above functions would be advantageous, especially when the ex- 
pansion is being used for sub-tabulation. For purposes other than sub-tabulation, 
values of the Chebyshev functions can be obtained by simple interpolation in the 


existing tables [6] and [8]. For convenience, values of cos are given to 


os 
4j+2 
twelve decimal nlaces in Table II. 
G. A. CHISNALL 


Department of Astronomy 
The University 
Manchester 13, England 


I am indebted to Professor Kopal for allowing me to read the relevant sections 
of his book prior to publication and for many constructive criticisms. This note 
has been written during the tenure of a scholarship from the N.R.D.C., to whom 
I am also indebted. 
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ee te ie, Sect. ‘A, v. 62, 1946, p. 204-210. 
. E. W. BROWN, Tables of the Motion of the Moon, Yale Univ. Press, v. 1, {pine 110. 
K. Camp, “The use of least squares in interpolation,” Actuarial Soc. of 'Am., Trans. -» V. 29, 
1928, p. 216-225. 
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5. L. J. Comrie, “Construction of tables by interpolation,” R Astr. Soc., Monthly 
—_ v. 88, 1928, p. 506-523. ai 
ry : MILLER, op. cit., p. 196. 
‘ ‘ ASSACHUSETTS INSTITUTE OF TECHNOLOGY, Mathematical Tables Symposium, Cambridge, 


. NBS Applied Mathematics Series No. 9, Tables of Chebyshev Polynomials S,(x) and C,(x), 
U. S Govt. Printing Office, Washington, 1952. 

C. W. CLensnaw & F. W. J. O_ver, “The use of economical polynomials in mathematical 
tables, ” Cambridge Phil. Soc., Proc., 51, 1955, p. 614-628. 

Note added in proof. Dr. L: Fox of the National Physical Laboratory, Teddington, Middlesex, 
England, has recently informed the — in a paper as yet unpublished, that in fact the Everett 
form of the modified Chebyshev-Bessel interpolation formula gives slightly better results (in the 
sense that the maximum possible residual error is smaller) than the toon & advocated here. 


Iterative Procedure for Evaluating a Transient 
ResponseT hrough its Power Series 


Introduction. We discuss here a particular method of evaluating a time 
function, such as a transient response, from its Laplace transform. We shall 
assume the Laplace transform of F(?), 


(1) P(S) = f ” e*F(t) dt 
0 
is given, and is of the form 


Ags’ + As’ +---+ A°y_2S + A°y—1 
SY + byS*—! + b.S%-? +---+ by_iS + by 





(2) P(S) = 


(The superscript 0 of the A’s does not indicate a power, but is used as a super- 
script for reasons which will appear presently.) 

We require in (2) that the denominator be of higher degree than the numer- 
ator, which is necessary if F(#) is a Laplace transform. We shall always assume 
also that we have F(S) written with the coefficient of S¥ in the denominator equal 
to unity, as in equation (2). 

The method described below is an iterative method for obtaining the quan- 
tities a; in MacLaurin series for, F(¢), which we write in the form 


(3 FP) = Day 
imo )=J- 


The method has the advantage that it is not necessary to know the roots of 
the numerator or denominator of F(S). The series (3) will converge for all 
positive values of ¢, although for larger values of ¢ a great many terms may be 
needed. This, however, is not too great a disadvantage for digital purposes, 
since the method for computing the a;,'s is an iterative one. It is necessary, how- 
ever, to have some idea of how many terms will be needed, and a criterion for this 
is also given below and discussed. The theoretical derivations are given in 
the Appendix. 

Method for Computing the a,’s. The method for computing the a,’s, the theory 
of which is discussed in the Appendix, is as follows: We start with the coefficients 


A, A’, nite A°y-1, bi, ba, es by of (2). 
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Then 
ag = Ao. 
We next evaluate the numbers Ao', A2', ---, A'y_; as follows: 
Aj =A; — abi 
A} =A; — ao: 
(4) 
A'y-2 = A®y-1 — aoby-1 
A'y.1 = — aoby. 
Then 
a; = Ajo. 


Then we evaluate 
Ag = A? —_ ab; 


A; = A? — abe 
(S) 


A’y_s = A'y-1 — aiby-1 


A’y_1 = — aiby 


and a: = Ao’, and so on. (Again the superscripts 1, 2, etc., are merely super- 
scripts and do not indicate powers of the A’s.) In contracted notation we have 


a; = Ao? 
(6) Aj! = Adar — ajdeys, OS RS N-—2 


when j = 0, 1, 2, 3, etc. 

Thus the numbers a; can be obtained and substituted in equation (3) to 
give F(t). 

Number of Terms. It is difficult to estimate the exact number of terms of the 
series (3) necessary to give a prescribed accuracy. However, it is shown in the 
Appendix that if we wish to have F(#) accurate to within +e with 0 < ¢ < hy, 
neglecting numerical round-off errors, it is sufficient to have nm > mo, where 
(n + 1) is the number of terms (counting the 0 order term), and 





N—4 
(7) my = e(1 + dy + In| 40 Slee tite 1], 
2be(e — 1) 
here A is the maximum of the numbers |Ao°|, |A:°|, ---, |A®°v-a|; 5 is the 
maximum of the numbers |6;|, |52|, ---, |dw|; and e is the base of the natural 


logarithms—‘“‘In” denotes the natural logarithm. 
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It should be borne in mind that the above criterion may give more terms 
than are actually necessary in a particular case. However, it has the advantage 
that it offers a method of determining the number of terms from the function 
F(S) alone, without necessarily having advance knowledge concerning the func- 
tion F(t). 

As an illustration we consider two simple cases where F(t) is easily evaluated 
by elementary methods. First consider 


1 
MS = S47 
Then we know 
F() = et. 


(Our iterative procedure will give a; = (—1)/.) 

In this case, A = 6 = N = 1. If we take t = 5 seconds, and « = 10-*, sub- 
stitution in equation (7) gives m) = 31.3, so we would use 33 terms (counting 
the 0 order term). However if we take = 17, our error is less than 5'*/18! which 
is less than 10-*. 5!7/17! is greater than 10-*. Thus the actual number that we 
need in this case is m = 17, or 18 terms, so that equation (7) gave us far too many. 

Next consider 


2S + 11 


PS) = Sy iis +10 





(N = 2,A =b = 11). 


Here again we can find directly that 
F(t) = e*+e-™ 


so in this case 


a; = (—1)i[10i + 1]. 


Substitution in (7), with the same ¢ and « as in the previous case, gives 
m = 171.6 so that we would use 173 terms according to our formula. The 
number of terms that is actually needed in this case is around 145. Thus in this 
case again our formula gives too many terms, but not by as large a ratio as in 
the previous case. 


APPENDIX 


The method under discussion is based on the expansion of F(S) around the 
point at infinity. Thus we substitute 1/Z for S and get 





(8) F(Z) = P( s) = Z{A0® + AYZ +---+ A*wiZ"} 


Zz 1+ b:Z +-+-+ byZ” 


Next we write 


(9) f(Z) = Zg(Z) 
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so 
(10) ez) - 
where 
(11) po(Z) = Ae’ + AYZ +---+ AwiZ*4 
and 
(12) Q(Z) = 1+ b:Z +---+ byZ*. 
Now g(Z) is analytic at the origin (Z = 0), and so it has a MacLaurin series: 
(13) g(Z) = ao + aiZ + a2oZ? + ---. 


We shall show presently that these coefficients a;, are the same as the a; in 
equation (3), and also that the series (3) converges for all values of f. 

To evaluate the a; in equation (13) we use the method of [1], which becomes 
somewhat simpler in this case since we are expanding around the origin. From 
(10) and (13) we have 


(14) bo(Z) = Q(Z) [ao + aiZ + a2Z* + ---). 
Then 
— 60(9) _ 4 
(15) a Q(0) Ao . 


Now, having ao, we note that 


bo(Z) — acQ(Z) 


has 0 as a root, so we can define another polynomial 


po(Z) — asQ(Z) 





(16) pi(Z) = z 
Then 
(17) pi(Z) = Ad + AYZ + AZ? +---+ AlyiiZ* 


and the coefficients A,! are given by equations (4). We then have 


(18) pi(Z) = Q(Z)[on + asZ + a2? + ---J, 
Thus we get 
— 2100) _ 4a 
(19) a Q(0) Ag . 


Then defining the polynomial 


p1(Z) — a0(Z) 


p2(Z) = Z 
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we get 
(20) 2(Z) = Ae + AvZ +---+ A*yiZ*" 


where the coefficients A;* are given by equations (5). Then a: = A¢’. If we 
continue in this manner we get the iterative relationships (6). 

Next we shall use the method of [2] to prove equation (3), since this will 
assist us in obtaining an error criterion. Suppose « is a small positive number. We 
want to show that in a given interval, 0 — t, we can choose m» such that when 
n> No, 


n ti 
lr = Laz <6 O<h bh. 
im ~=s J? 





First let 7; be a positive number such that all the roots of Q(Z) are greater than 
r; in absolute value. Then g(Z) is analytic for |Z| < 11, and so also is 





_ pi(Z) 
(21) h(Z) = 0(Z) 
From equation (18) we see that 
(22) h(Z) =¥ az. 


j=1 


The series (22) will converge uniformly for |Z|<1r1. Thus if we let 


(23) RZ) = E az 


j=ntl 


then for any positive number ¢; we can find mo(¢:) such that |R,(Z)| < « when 
|Z| <r: and m > mo(e:). 
€, can be any positive number and we will choose 


ee 


(24) =e ", 


v 


1 
Since the point — + 10 is to the right of all the singularities of F(S), we have 
1 


by the inversion theorem 


i... 
=—+4 a 


(25) 2niF(t) = f i P(S)e"as = f a en | = | as. 





a, eee 
=e i=0 
Lal 


The last step follows from equations (8), (9), and (13) substituting 1/S for Z. 
Next we consider the function 


(26) F,,(t) = Sen: 
imo |= J 
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1 





























The Laplace transform of #/7! is om’ (j = 0,1, 2, ---). 
Thus by the inversion theorem 
27) 2riF.) = Lai | f” <as| 
( n ia -] a. Siti ” 
Then for m > mo(e:), we have 
(28) |F() — F.(| 
ae Sane ae! 
aj 
a3 oom et i ‘. aagnimenn 
on - i 2 on dS Xai : Sit as| 
1 Ra i : t 
1 ate R(3) + ixt ates 1 
tl ‘ oki tere : ‘- 
al. “Ss SI (aay | 
Jaa 7 + 
rT 
4 tty 
ere d " 
Tr 0 T 
Fein (55 
T1 


if 0 < ¢ < t. The numbers, #) and « may have any positive values. Once they are 
chosen we choose ¢; by equation (24), and then find mo(e:). The actual finding of 
a number mo, after having found ¢; is the task which we will consider next. How- 
ever we have proved that the series of equation (3) converges to F(é) for all 
positive ¢. Also we see that the number m which makes | F(t) — F,(#)| < «¢ for 
t = to, also does so for all values of ¢ between 0 and ¢. “ 

Next we ask the question: Having given ¢, from equation (24), how do we 
actually find a number mp such that |R,(Z)| < «: whén |Z|< ri and m > m? 
We must do this in order to have a formula for obtaining a sufficient number of 
terms in the series (3). Also we have not yet chosen 11, except to say that it is less 
than the modulus of each of the roots of Q(Z). 

First of all we define 


1 
2 = — 
(29) ~~ oo 
where b is the maximum of the numbers |}:|, |d2|, ---, |dw|. Then for |Z| < ro, 





(30) jo@lz1—s[454(545) +--+(G5)'] 
-1-( b yt) 1 
1 


1=Eeins) ~ +5)" 
1+ 





o 





Thus 


(31) 


wher 


(32) 


(33) 


wher 
we h 


(34) 
and 


(35) 


if |Z 


(36) 
The | 


(37) 


(38) 
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Thus all the roots of Q(Z) are outside the circle |Z| = 7». Now we choose 
6 
(31) "n= Oro = 1+5 ® 
where @ is as yet undetermined, except that 
0<@é@<1. 
Then for |Z| < ri, we have 
(32) R.(Z) = YL aZ* = h(Z) — Da Zh. 
jenti j=l 
Also 
_ h&-2(0) 
(33) %= Gr! 
where h®(Z) is the k® derivative of k(Z). From Cauchy’s integral formula, 
we have 
1 £ bw) 
(34) h(Z) = ~=t ¥ - zw 
Se and 
k! h(w) 
35 h®(Z) = — ——ae Gay 
are 5) ) 2mi Sc, (w — Z)*# 
g of ; , : 
aes, if |Z| < r1, where Cp is the circle | w| = ro. 
- all Thus, for |Z| < 71, 
for 1 h n Zi-l h 
(36) R,(Z) = $ 2O ww - o = ¢ Maw 
2riJo,w—-Z ju 282 Jo, wi 
> we 
Me? The second term 
er of 
; * Zr h 
™ len DS g ©) tn 
jai 271 Jo, wi 
(-2) 
1 h(w) 2 (Z\>° 1 h(w) w 
=— —— = dw = — —— ——— dw 
2xt Co w j=1 w 2x1 Co w (1 aes 2) 
< To, bad 
ae a. I dw — ae wa (2) dw 
2niJo,w —Z 2xi Jc, (w — Z)\w ’ 
Hence 
nv 1 h(w) (2 y 
b)* owes = 
(38) R,(Z) il Ja @- Sj \e dw 
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Thus, since |Z| < ri = Oro, 
1 M Me 
39 R,(Z)| << — ———— #2er, = ——— 
(39) IRaZ)1S a aw 
where M is the maximum of the function |#(w)| on the circle Co. 
On the circle Co, we have 


| Pi(w) | 


(40) (A) = TG 





< (1 + 5) |pi(w)|, 


using (30). 

The coefficients A;' of ~:(Z) are given by equations (4), and since ao = Ao’, 
we see from these equations that 
(41) |Ai|< A(i+5), k=0,1,---,(N — 1), 


where A is the maximum of the numbers | Ao°|, |Ai°|, ---, | A®w_1]. 
Thus on the circle Co, 





42) lows aa +o [1+rt+( 5) 4--4+(4)"] 


N 
‘- (5) 
1+) 


Oe cineeee 


1+6 





= A(i + d) 
Therefore, using (40), 
( 1 y’ » 
a on _ A(t + bE (1 + 6) — 17 


papi BN 6 
1+6 


M<A(i + 0)¥* 





1 


Thus, returning to condition (39) we see that for |Z| < rn, 


A(i + b°C(1 + 5)¥ — 1)" 


(43) |Ra(Z)|< ha 





Thus, we wish to have 


A(1 + 6) (1 + 6) — 1] _ 
b(1 — 6) ™ 





This will be true if 





1\"_ A(i + 5)*{(1 + 5)¥ — 1] 
(5) “ bes(1 — 6) 
that is, 





A(i + b)*[(1 + b)¥ — 1] 
in| be:(1 — 8) | 


“() 





n> 





Thu 


(44) 


(45) 


(47) 


whi 


Hor 
to € 
mat 


A’, 
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Thus we have found a number mo(e1) which has the required property, namely 





A(1 + 5)[(1 +5)" — 1] 
“ [ ba(1 — 6) | 


=) 


In order to relate ¢, to «, the desired maximum error in F(#), and the interval 
(0 to to), we use equations (24) and (31), and get 





(44) ny = 


+b) 
- [ 40 +800 + ee ’ | 
2b(1 — Oe 


(3) 


A(1+5)[(1 + 5)¥ — 199 
a+o, ™ [ 2b(1 — B)e | 


“a(’) aC 


We still have not chosen 6, except to say that it is a positive number less than 
unity. We somewhat arbitrarily choose @ so as to make the f) term as small as 
possible, i.e., the value which makes @ In (1/6) maximum. This value is 








(45) 








_ 





(46) thea 
Thus 
eae es 


which is the same as equation (7). 

This paper was written while the author was employed by the Minneapolis- 
Honeywell Regulator Company of Minneapolis, Minnesota. The author wishes 
to express his appreciation to Mr. Eric Clamons of Minneapolis-Honeywell, for 
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The Numerical Integration of Two-Point 
Boundary Value Problems 


1. Introduction. A general system of linear, ordinary, differential equations 
may be written, in matrix notation 


(1) Y =A(t)Y + F(t), 


where the dot denotes differentiation with respect to the independent variable, 
t; Y is the column matrix of unknowns, 


¥1 


Vn ’ 


A(t) is a square matrix, and F(é), the forcing functions, is another column matrix, 


fi 


fn 

Suppose that r values y,(0) --- y-(0) are given and at ¢ = 7,(m — r) of the 
y's are given, viz.: y;,,& = 1, ---, (w — r). In order to solve (1) using high speed 
computing machinery, whether analogue or digital, it is essential to have n condi- 
tions at one end of the interval 0 — T. One possibility is to carry the constants 
¥r41(0) --- yn (0) along as unknowns during the integration, and then to determine 
them when the far end is reached. Another possibility consists of guessing 
¥r41(0) --- yn(O) and then adjusting them. These techniques have been described 
by Tifford [1]. The first of these methods is awkward especially if a large number 
of unknowns must be carried; the second requires more integrations than the 
methods described below. 

This paper is concerned with methods for converting, a two-point boundary 
value problem into an initial value problem. For linear systems the conversion is 
accomplished with a finite number of integrations. In section 2, where nonlinear 
systems are considered, an iterative procedure is required. 


2. Linear Systems—Method of Adjoint Equations. The adjoint system was 
developed by Bliss [2]; more recent developments are described by Wright [3]. 
The system adjoint to (1) is, by definition, 


(2) -X = A(x, 
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where the bar denotes the transposed matrix and X is a column matrix 


xX) 
Xn 
The x,’s and y,’s are related by 
(3) © (egy) = aif 
dt Wi = Xijsi, 


where the repeated suffix implies summation from 1 to . A proof of (3) follows: 
Consider 


d d : 
7 =F (XY) = XY + XY. 
Substituting for the derivatives from (1) and (2) there is obtained 


¥ (eg) = RAV +H) - (ANY 
= YAY+XF-— (XA)Y. 


The first and last terms cancel identically and the result (3) follows. 
Integrating (3) over the interval [0,77], it is found that 


(4) 8 xi(T)y.(T) — x4(0)9-(0) = f ” walt) Fulda. 


This fundamental result is the one-dimensional form of Green’s theorem. 

The aim is now to calculate y,,:(0) --- y,(0) using the given conditions. The 
method described in this section will utilize the adjoint system (2), whereas in the 
next section the original equations (1), only, will be used. 

The adjoint system (2) must be integrated (» — r) times from T to 0 subject 
to certain starting values. These are chosen such that the x;(7)’s are all zero except 
one. For the m** integration this exception is the x; which appears as the coeffi- 
cient of y;,(T) in (4), and which may be designated by x;,(7). For simplicity 
x;,(T) is chosen to be unity. 

Let »x;(#) be the value of x; deduced from the m* integration of the adjoint 
system. The fundamental result (4) holds for each integration of the adjoint sys- 
tem, and this is expressed by rewriting it 


6) wxi(T)yi(T) — axe(0)9s(0) = f ” ealt) filE)dt. 


But since all x;(7)’s were chosen zero except x;,(7), the first summation reduces 
to a single term and (5) becomes 


6) wx:(0)9:(0) — y4,(T) = L wxi(t) felt. 
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This is a system of (m — r) linear algebraic equations for the (n — r) unknowns 
y(0),t=r+i---n. 

With these values determined (1) can now be integrated from 0 to T, as an 
initial value problem. The total number of integrations consists of (n — r) of the 
adjoint system plus one of the original given system. 


3. Linear Systems—Method of Complementary Functions. In this case the 
homogeneous equation, derived from (1), 


(7) u = A(t)u 
is integrated from 0 to T (nm — r) times. The m** time the initial conditions are 


nts(0) = 0 ix¥r+m 
1 t=r+m, 


ll 


where ,,u;(¢) is the result of the m** integration of u;(f). 
Next the non-homogeneous equation 


(8) V =A(t)V + Fit) 


is integrated once, from 0 to 7, subject to 


S sace 
r+i--- am. 


v:(0) = »:(0) 7 
= 0 1 


The general solution of (1), subject to the given conditions can be written 
(9) yi(t) = Cy pus(t) + 0, (4), 

where the C,’s are constants determined from the (m — r) equations 

(10) yu(T) = Cp p(T) +24,(T) k= 1--- (mn —71). 


In this equation the left-hand side is given and the quantities u;(7) and »;(T) 
have been determined from the integrations of (7) and (8). By letting ¢ = 0 in 
(9) it is seen that Cp, = y-47(0), Pp =1---n—r. 

The method of complementary functions takes less time than the method of 
adjoint equations for a single problem, but it is not so suitable when solutions are 
required for a large number of different forcing functions, F, or for a large number 
of different boundary conditions. 


4. Nonlinear Systems. A general system of nonlinear, ordinary, differential 
equations may be written 


(11) Vit) = gilyr +++ Wnt) T= 1--- 0, 
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where it is assumed that the functions g; are once differentiable with respect to 
all the y,’s. Once again the boundary conditions are specified as y:(0) --- y,(0) 
and y;,(7T), k =1--- (n — 1). The procedure is to estimate y,,:(0) --- y,(0) 
which estimates are denoted by y*,4:(0) --- y,*(0). Equations (11) can now be 
integrated as an initial value problem, and the solutions are denoted by y,*(é). 
In general, it will be found that the computed values of y,;,(7),& = 1 --- (n — r) 
namely, y*;,(7), differ from the given values of y;,(7). To obtain the correct 
solution of the system (11) the difference between y*;,(7) and the given y;,(T) 
must be made as small as possible. 
Define éy;(t) by 


(12) dy:(t) = 9: —y*@) t= 1--- 0. 


Substituting y;(#) in (11) there results, to a first approximation, 


ogi 
(13) a3.(0) = (2) by 


The equations adjoint to (13) are 


(14) —X; ~(2) =. 


Since non-homogeneous terms are absent from (13) the Green’s theorem for (13) 
and (14) reduces to 


(15) xi(T)dyi(T) — x;(0)8y.(0) = 0. 


The quantities required are dy,(0), i =r+1---m. To obtain (nm — r) 
equations, for these (m — r) unknowns, (14) is integrated from T to 0 (m — r) 
times; the result of the m** integration is denoted by ,,x;(#). Equation (15) holds 
for each integration of the adjoint system, and this may be expressed by 
rewriting it 


(16) mX<i(T)dyi(T) — mxi(O)dy,(0) =O m=1--- (n—71). 


The m** time the adjoint equations are integrated the starting values are that 
all x;(T)’s are zero except the coefficient of 5y;,(T) which is chosen to be unity. 
The first summation in (16) now reduces to a single term; furthermore, éy;(0), 
| = 1---,r are identically zero. Hence (16) reduces to 


(17) 5yin(T) — mXryi(0)5y-4:4(0) =O m=1--- (n—7), 


where the summation on 7 extends from 1 to (m — r). Substituting these values 
into (12) the improved estimate of y,;(0), i = 7 + 1--- m are obtained. In the 
case of nonlinear systems (11) is integrated again, using the improved estimates, 
and the process is repeated. In the special case where the g;,’s are linear functions 
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of y; (13) is exact and the new estimates of the initial values are precise after 
only one application of the method. This method, when applied to linear systems, 
is similar to the second method described by Tifford [1]. 


5. Remarks. The problem was posed such that the y,’s were given at either 
end of the interval 0 — 7. The methods are easily extended, if instead,  — p 
linear relations are prescribed between the y;(0)’s and » between the y;(7)’s. 
The methods may also be generalized to convert three or more point boundary 
value problems into an initial value problem. 

An alternative method, suitable for nonlinear systems, has been presented by 
Lance and Deland [4]. This method is also one of iteration but it involves numeri- 
cal differentiation and for this reason is less accurate per iteration than the 
method presented in section 4. 

If large scale computing machines are not available then relaxation methods 
can be used to solve systems of differential equations. Allen [5] states that if 
relaxation methods are to be used then an equal number of boundary conditions 
must be specified at each end of the interval, and for this reason the number of 
equations, ”, must be even. Allen presents a method for converting problems which 
are not stated in this form into this form. The authors considered using the 
methods given above as an alternative to Allen’s. Although, it is possible to do 
this, Reid [6] has shown that the conversion is unnecessary and that relaxation can 
be used directly regardless of how the boundary conditions are distributed. Reid 
also shows that the number of equations need not be even. 


T. R. GoopMANn 
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High Accuracy Quadrature Formulas from 
Divided Differences with Repeated 
Arguments 


In the numerical solution of differential equations, one uses quadrature 
formulas to relate the unknown functions to their derivatives. While the functions 
enter only at the beginning and end of the range of integration, all intermediate 
values of the derivatives may enter. For checking purposes at least, it would be 
advantageous to increase the accuracy of the quadrature formulas by permitting 
weights to be assigned to intermediate values of the functions, as well as their 
derivatives. An interesting method of obtaining such formulas is given below. 

Given a function f(x) and m + 1 distinct values of the argument xo, x1, - - -,-<a, 
the mth order divided difference of f(x) with respect to these values of x may be 
defined as 





”. f(x») 
1 ,X1,°**,X_) = : 
(1) f (xo, x Xn) ~ IL. 
where 
(2) II, = II (&»— =); 
for example, - 
f (xo) Ff (x1) f (x2) 


(3) f(xox1,%2) = 





(xo — x1) (x0 — x2) | (x1 — x0) (x1 — 2) | (x2 — 0) (X2 — 21) 


An upper limit can be obtained for the value of a divided difference by means 
of the following theorem, see [1]. If f(x) and its first » — 1 derivatives are finite 
and continuous and the nth derivative f‘” (x) exists, then 


# ©) 


n! 





(4) Ff (xo, *1, +++, Xa) = 


where ¢ lies in the range between the smallest and largest of the arguments 
Xo tO Xp. 

For a divided difference involving a repeated argument one takes the limit as 
two initially distinct arguments approach one another; thus 


Ff (%o’) fo) | _ @ 
| —- & + xo | & dx, 1 **)- 





(5) f(xo, x0) = Lim f(x0’, x0) = Lim 


In general, as shown in [2], 
ko +1 ki+1 ke +1 
—_—_—S eee 


-__, 
(6) F (X0,-+ ++, Xo X1, *** 5 Xty Muy °° *y Xn) 


1 the thi t-:-+he 
~ Bolkal---Bal Oargday™ «~~ Oxrgbe 100 *H °° Fo) 
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and, for the special case k; = k, for all i, 


k+1 k+i1 k+1 
oN oe —_ 
(7) F, = f (xo, **% 5 XO, Xi, °° *y Kip **%y Kn, 72+, Xn) 


Wasct OF | fl») 
ins apa OS | z a 
where JT, is given by (2). 


From (7) 


-—!_ | 7% | f@e) 
Oi: Bao aa = [1] Il, 
i 


~ Gp E ea fee LH aa i 


ixp 





and since in (8) 


(9) wt, ee eT) for 1p, 
one has 

k+i1 k+1 
(10) Fy = f (Xo, ++, %0, °° *; Say ***;%0) = 
Therefore by Leibnitz’s Theorem 


k+1 k+1 
Sl ameaneetnenmmeet n 


——~ k 
(11) f (x0, ++, X0, °°, %ay °**) Xa) = DON, Ff (x,), 


p=0 s=0 - 








where f‘*)(x,) represents the sth derivative of f(x,) and 


i 
s!(k — s)! dx,*-* JJT,*4 





(12) N,' = 


Special Case k = 1 
From (12) 


ti) 
(13) NS = >I? = - 211, 
Pp 


where 


(14) iis’ = - II, 2 
‘33 








Xp — Xj 


For equally spaced arguments x; = x9 + ih, where i = 0,1, ---,m, one has 
from (2) 


(15) II, = (—1)""*p!(n — p) th. 





(1 


(1 
th 


a nh 2h hCUe [06 oe See 
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Likewise, from (14) 








(16) II’ = (—1)*-*p!(n — p) ICS, — S,_p]h™", 
where Sp = 0 and 
(17) 5S, = Es, for r +0; 
therefore, from (13) 
No = —2LS, < S.-»] 
> ~ PALp ln — p) IF 
(18) 
N,' = : . 
” — ‘W[p!(n — p)!P 


Multiplying (11) through by (m!)? h?*+" and making use of (4), one has 


ntl 


n n h? 
(19) LD Ansf (Xp) = h LD Bagf (xp) — D 
p=0 p=0 





fer (é), 


where xo < ~ < x», and one sets 


(20) Anp 2LS>, - S.-»][C,"*F 
(21) Bay = (C,"F 
and 

_ (Qn +1)! 
(22) D, = ~~? 

In the above equations 
n! 

i Co" = Sin — DD! 


are the binomial coefficients. 

The first two quadrature formulas corresponding to m = 1 and m = 2 are the 
trapezoidal rule and Simpson’s rule, respectively, applied to the derivative of 
f(x). The formulas for higher value of m have the usually undesirable feature of 
involving the value of the integrated function f(x) at intermediate values of x 
between xo and x,. They should, however, due to their very small remainder 
terms, prove useful in checking the accuracy of integration in the numerical 
solution of differential equations. Moreover, these quadrature formulas should 
be advantageous in the solution of two-point boundary value problems, since 
these must be solved by a relaxation or by an iteration method. 
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Tables of the coefficients in (19) are given below: 


-e SS SS SS eS OC 


ork Se 


25 
36 


— 


36 
100 
225 


TABLE 1 
Coefficients Any 
2 3 
3 
5 i 
0 2 
_ 200 200 
3 3 
525 
“4 
TABLE 2 
Coefficients Buy 
3 4 
1 
16 1 
100 25 
400 225 


36 


1 


137 
30 

462 
5 





49 
10 


140 
630 
2772 
12012 
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On a Device for Computing the e,.(S,) 





Transformation 
1. Introduction. The transformation 

a Sati 5 ins } ae 
AS, ASa41 +--+ ASn4m 
AS, °° AS. se, 

(1) em (Sn) — +m—1 + 2m—1 
1 1 ---1 
AS, ASa41 —— ASaim 
AS n4-m—1 7 ASn42m—1 








has been developed by Schmidt [1] and Shanks [2]. It is shown by both authors 
that 

éu(Sa) =a 
if 


(2) Si =a + = b-e,". 
r=1 


The former writer shows that iterates of the form (2) occur in the Gauss-Seidel 
method for the solution of linear equations, and derives the condition 


AS,, ASn41 +++ ASntm 
ASn+1 ASn+2 peek: ASn4m+1 
(3) es Ss - «4 6 eee 
ASaim ASa4m+1 °° ASa+2m 
for the applicability of the transformation. 

A consideration which militates against the use of the transformation is that 
the vast amount of labor spent in evaluating the determinants in (1) serves 
only to produce one transformed result. The computation of further results pro- 
ceeds quite independently with a similar expenditure of effort, so that use of (1) 
in its present form as a sequence to sequence transformation in which members 
of the transformed sequences ¢,,(S,)(m = fi, i + 1,---,m =1,2,---) may be 
examined for any tendency to approach a limit, involves a prohibitive expenditure 
of labor. It is proposed to show that the transformation (1) may be effected by a 
simple algorithm, in which transformed results for m = %, fi + 1, ---,m = 1,2, --- 
are progressively available for comparison. 


2. Theorem. Jf 
€om (Sn) = €m(Sn) 


and 
i 
€2m+1(Sn) = en(ASu) 
then 
(4) €041(Sn) ad €.—1(Sn41) + : s=1,2,---, 





€s(Sn41) = €.(Sx) F 
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provided that none of the quantities €o,(S,) becomes infinite. 
Proof. The above result is easily verified when s = 1, for: 





1 
= S, —_____ 
€:(Sn4i) — €1(Sa) ar | 1 


Atus 45, 
Sn41(AS, — ASaz1) + ASa41°- AS, 
AS, — ASn+1 
SrASn41 br 41° AS, 
ASn+1 — AS, 
€2(S,). 


€0 (Sa+ 1) + 











For larger values of s, consider first the case s = 2k. It is then to be proved that 











1 1 ++ 1 1 1 os | 
A?S, A*Sa41 paola A*Sas Sees Sus os A*Suss 
5) A*Sa42-1 eee A*Sn+2e-1 A wane eee A°Sessn-s 
AS, ASa41 a ASn+k ASn+1 ASn+2 = ASa+t 
ASasz1 ASaz2  -°* ASapegi ASa+2 ASais °°: nae 
ASn+k ASa+e41 pied a ASn+2% ASa+% ASn4e41 aia te Mins 
1 
Sa+1 Sn+2 7 Sapey1 S, : Sa41 ses Satt 
ASn+1 ASa+2 tts ASasae1 AS, ASn41 sic ASn+t 
ASaiz ASniey1 -** Swan ASnze-1  ASaze +++ ASap 2-1 
1 1 ae | 1 1 oes | 
ASa+1 ASni2 -** ASesaes AS, ASn+1°** ASees 
ASnie ASn+e41 °° a Sea ASa+e-1 cee es 














Rearranging columns and rows, the left hand side of equation (5) may be written 

















1 1 ess 1 1 i pan | 

OSes. A*Sai2 st A*Sa+t A’S,, aSuse wSus a A*Sn+k 
MSaa A*Sn+k+1 eee A*Sn+2%-1 A*Sn+4e-1 BScaa1 coe A*Sa+24-2 
ASa+1 ASn+2 sata ASa+k AS, ASn+1 ASn+2 — ASa+t ; 
neal Sue ik a A*Sa+k A’*S, aes oS: > A*Saue 
A*Seus A*Seanys -=* A*Sason-1 A*Saye-1 A Seat A’Seas » A*Sa428-2 

















at 
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and by the Schweinsian expansion of the quotient of two determinants [3] this 
expression is equivalent to 























1 1 — 1 1 A*Sni1 A*Sa+2 ” A*Sast 
ASa+1 ASa+2 = ASa+x AS, A*Sa+2 A*Sa+3 7 AS44241 
A*Suss APSugs - A°Saiz A*Sp i diaBe taaietd cf  R 
(6) Sorat A°Sese ae A*Sn+2%-2 A*Sn+e—2 A*Sa+t A*Sn4e41 > ait 5 A*Sa42%-1 
ASn+1 ASa+2 “— ASn+ AS, ASn41 ASa+2 ee ASast 
wuss baeen ad A*Sa+k A’S, au A Suse “2Se00 
A*Svss BSesser +++ A*Sasoe-1 A*Sn4e-1 A*Seas-1 A*Seas - -A Seis 
1 1 -+-1 1 1 ---1 
a aoe so atoes ASn+1 ASn+2 arty ASn+e+1 
ASoss-1 ASwus “ae pret ASa+t ASn+e41 see ASn+i2% 
AS, ASa+1 > ASa+e ASaz1 ASa+2 > ASa+e 
ates hres » ASn+e¢1 ASa+2 Shes > AEoens 
ASwsa ‘ASeeass — ASn+2% ASa+e bea go han 
The right hand side of equation (5) may be written 
1 1 --- 1 1 1 -1 
(- 1)* AS, ASa+1 “ita ASn+t aSirs BS i ASn+e+1 
(7) ASnze-1 ASaye +++ ASny2n-1 ASuas pe +++ ASnise 
aes ane > ASa+e+1 ASa+1 *** ae AS, 
ASwss pe wh Bp ASn+2% ASa+t ee jets b ASa4e-1 
Sati Sage > Sn4ept 1 oe | 1 
ASn+1 Pe? ates AS, ASni1 Sore "Aiea 
~ | ASage ++ ASexse-1 ASn+e-1 ASn+e _ aw °° ee 
Sagi *** Sate Se 1 1 ---1 
and by using an extensional identity derived from 
a, de 
lai] |b2|—|d:| | ae| 


bi be 
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(7) is equivalent to 





(—1)* AS, ASa41 ? ASn+t ASa+1 ASa+2 oor ASn+e41 
ASn+e+1 ASn+e he his ASn+2%-1 ASn+e ASa+e+1 Fg om ASn+2% 

ASa+1 ASa+2 nee ASa+t ASn+e+1 AS, 

. ° = % ‘ ° . ° ° ° ° ° ° ASn+1 ASna+2 -++ ASaae 
ASn+2 ASais +++ ASnyey1 

ASaze ASnreei +++ ASnpor-1 ASnzee ASa4e-1 eae pepe, eye) * 

Sati Sate +++ Suse Sattti Ss ASnse ASnzep1 ++ ASny2n-1 

1 1 ae | 1 1 








which is equivalent to expression (6). Thus the left and right hand sides of (4) 
are equal when s = 2k. That they are also equal when s = 2k+1 is demon- 
strated in a similar fashion. 


3. Algorithm. It can now be seen that the quantities S,, «,(S,);” = #,% + 1, 
+++; = 1,2, ---; may be arranged as in a tabie of S, and its differences, in which 
¢,(S,) takes the place of A*S,. This table is composed of groups of four quantities 
€(Sa), €2(Sn41), €e—1(Sn41), €241(S,), arranged in the form 


€s (Sn) 
€s—1 (Sn+ 1) €s+1 (S,). 
€s (Sn41) 
The quantity ¢,,:(S,) is formed by the relation 


1 
€2(Sn+1) _ €.(S,) : 





€24-1(Sn) = €s—1(Sn41) ss 


and it may be remarked that this simple calculation is easily programmed for a 
digital computer, all the entries in the second and subsequent columns being 
obtained from this programme by trivial alterations in storage parameters. 

In this table the even order columns €2,(5,) display the transformed results 
e,(S,). There is no need for the explicit evaluation of the determinant in equation 
(3), as it occurs in the table as the denominator of ¢2n+41(S,) ; if equation (3) holds 
when n = fi, # + 1, ---, then €am(Sa) = €2m(Sazi1) = «°°. 


4. A numerical illustration. A numerical illustration is provided by the trans- 
formation of the sequence produced by the linearly convergent iterative scheme 


S41 ~ aa 4(S,? + 2) 


which may be used to derive the smaller zero of the Laguerre polynomial 
L2(x) = x? — 4x + 2, ie., x = 0.58578 64375. 
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It will be noted from Table I (and this is a consistent feature of the transforma- 
tion when meeting with numerical success), that 


1 
€2e(Sm+1) — €2s(Sm) 





> | €2e—-1(Sn+1) | 
and hence 
1 


1 1 
€2s(Smi2) — €2e(Sm+1) * €2e(Sm+1) — €22(Sn) 


Ceve(Sm42) — €26(Sm+1) ILere(Sm+1) — €2e(Sm)] 
[eee(Sm42) — 2eee(Sm1) + €2e(Sm) ] 





€2042(Sm) = + €2e(Sm41) 








= €25(Smii1) — 


It is thus an approximate rule that there are as many significant figures in the 
difference €2.42(Sm) — €2e(Sm41) as there are in the second difference ¢€2,(Sn+2) 
— 2e€ee(Smsi) + €2e(Sn). From this it is seen that loss of figures due to cancella- 
tion may occur in the transformation. This possible loss of figures in the Aitken 
transformation ¢,(S,) has been remarked upon by Olver [4] and Shanks [5]. 
Table I also displays terms in the transformed sequence ¢;*(S,), and it will 


TaBLeE I 
n S, € €2 €3 & e:? 
0 0.0 
2.0 
1 0.5 0.57142 85714 
16.0 89.11111 08 
2 0.5625 0.58510 63830 0.58574 34871 0.58576 19630 
60.2352 941 1658.713 
3 0.57910 15625 0.58573 19781 0.58578 57313 0.58578 61706 
211.05540 20262.3 
4 0.58383 96549 0.58578 18504 0.58578 64208 0.58578 64316 
725.9366 239060. 
5 0.58521 71856 0.58578 60462 
2483.836 


6 0.58561 97886 


be noted that comparison of the results given by the transformations e2(5,) and 
€:°(S,) in this case reacts slightly in favor of the latter. The difference between 
the two results however is not greatly significant, and there seems to be some 
point in using the transformation ¢,,(S,) consistently (in place of m successive 
applications of Aitken’s transformation) just in case S, might very accurately 
be represented by some equation of the form (2). There is no reason to believe 
that this is true for the example chosen, but results for a favorable numerical 
example are given in Schmidt's paper. Further examples in which success for 
ém(S,) might be anticipated from analytical knowledge of S, are given by Shanks. 


5. An application. Having shown that the various transformations e¢,,(5S,) 
may be effected without much effort, the primary purpose of this paper has been 
achieved. Functions of the form (2) are of such widespread occurrence in the 
theory of Mathematical Physics, that an exhaustive account of the problems in 
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which the transformation e¢,,(S,) might prove useful (besides being irrelevant to 
the main thesis of this paper which is primarily one of computational expediency) 
would be excessively discursive and, in view of possible developments, necessarily 
incomplete. By way of a footnote however it may be pointed out that the trans- 
formation may be applied to establish a criterion for the fitting of certain types 
of statistical data. The determination of the constants a;, b; when fitting tabular 
data, which is given at regularly spaced intervals in ¢ of magnitude w, to a function 
of the form 


rh 
(8) o(t) = X ber, 


when m is given, is already well known [6]. To establish the suitability of (8) 
as an adequate representation of the data, and further to determine 4, it is 
merely necessary to apply the transformation e,,(S,) to the data, writing 
Snir = O(t + rw), when e,(S,) = 0. (Various refinements, such as partitioning 
the data in groups of » values, and so on, are obvious and need not be discussed 
here.) This was one of the first problems upon which the transformation e,,(S,) [5] 
was used. 

It is the author’s hope that by demonstrating the ease with which the various 
transformations may be effected, their field of application might be widened, and 
deeper insight thereby obtained into the problems for whose solution the trans- 
formations have been used. 


P. WYNN 
Scientific Computing Service 
23, Bedford Square 
London, England 


J. R. Scumunt, “On the numerical solution of linear simultaneous equations by an iterative 
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pie. (2) and 49, formula (9). 

4, F. W. J. uve, “The evaluation of zeros of high degree polynomials,” Roy. Soc. Phil., 
Trans., v. 244, 1952, p. 385-415. 

5. D. SHanxs, “An analogy between transients and mathematical sequences and some non- 
linear sequence-to-sequence —— Ow —_— by it,” Part I. Naval Ordnance Laboratory 
Memento 9994, White Oak, Md., 

E. T. Warrraker & G.  iermdeiieg “The Calculus of Observations, Blackie and Son Ltd., 
4th ed, Ch. XV, 1944, p. 369. 


TECHNICAL NOTES AND SHORT PAPERS 


Note on the Computation of Certain Highly Oscillatory Integrals 


The purpose of this note is to draw attention to the possible use of the Faltung 
theorem for Fourier transforms as an aid to the computation of highly oscillatory 
integrals. 

If F is the Fourier cosine transform of f, defined by 


F(u) = Ef s0 cos ut dt 
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and G is the Fourier cosine transform of g, then 


[soe cos ut dt = if F(y)[G(u + y) + G(|u — y|) Jey. 


Numerical quadrature of the integral on the left will be laborious if u is large. 
But, if the transforms F, G are known, the integral on the right may prove much 
more tractable. The essential reason for this is that the large parameter u now 
appears in the argument of G as an additive constant and not as a multiplicative 
constant as in the integral on the left. 

As an example, consider the integral 


r=f (@ — 1)-* e-** cos bt dt 


for large values of 6. Here 





fd = eet F(u) = i e-w' lhe 


4 (2a) 


g(t) = (@—1)? for t>1 Gu) =- Ne Yo(u) 


=0 for ti<i1 
so that 


1 i f° 
git -iNES e-0L ¥o(b + y) + Yo(|b — y|) dy. 


This form is, for a wide range of values of a and 5, very suitable for computation 
using the simple trapezoidal rule [1,2]. The method was used very successfully 
in preparing a table of an integral, closely related to the one above, which occurs 
in the theory of the wave resistance of ships [UMT 174, MTAC, v. 7, 1953, 
p. 248). 


This note is published with the permission of the Director of the National 
Physical Laboratory. 


E. T. Goopwin 
National Physical Laboratory 
Teddington, Middlesex 
England 
1. E. T. Goopwin, “The evaluation of integrals of the form /_..” f(x)e~“dx,” Camb. Phil. 
Soc., Proc., v. 45, 1949, p. 241-245. 


2. Henry E. Ferris, ‘Numerical calculation of certain definite integrals by Poisson’s Summa- 
tion Formula,” MTAC, v. 9, 1955, p. 85-91. 


On a Cubically Convergent Process for Determining the 
Zeros of Certain Functions 


An iterative procedure of the form x»: = f(x») for the determination of the 
roots of the equation ¢(x) = 0 may be classified as linearly convergent, quad- 
ratically convergent, and so on, by using a notation due to Hartree [1]: If x = X 
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is a solution of ¢(x) = 0, x, = X + 9, (so that 9, is the error in the determina- 
tion of the root at each stage) and 7, is sufficiently small, then the convergence of 
the procedure x,4: = f(x.) is said to be of the k-th order if 9241 = 0(n.*). This 
means in effect that when using a quadratically convergent procedure we double 
the number of known figures in the solution at each stage, when using a cubically 
convergent procedure we treble this number of figures, and so on. An example of a 
quadratically convergent procedure is provided by the Newton-Raphson process 
o(%o) (where primes denote differentiation with respect to x). It 


$' (Xn) " 
may be shown that for this process lim Tet = ie) so that the Newton- 
oe Ga 2¢ (Xx ) 
Raphson process is indeed quadratically convergent. 

Various cubically convergent processes for extracting the roots of an equation 
are known, but they are not generally used because they involve the determina- 
tion of the second derivative ¢’’(x,) and the labor in doing so is normally con- 
sidered to offset the advantage of gaining the extra figures. It is the purpose of 
this note to point out that if the function ¢(x) satisfies a differential equation 
of the form 


Xn+1 = Xa 





(1) po" +46’ +16 =s 


(where ~, g, r, and s are functions of x), the determination of $’’(x,) is a trivial 
matter, once ¢(x.) and ¢’(x,) have been calculated. The number of known figures 
in X may thus be trebled at each stage with scarcely more effort than that 
necessary to effect the Newton-Raphson process. This remark applies equally 
well of course if ¢(x) satisfies a nonlinear second order differential equation, but 
the form of the latter permits of greater variety and hence the subsequent an- 
alysis is more complex. 

Perhaps the simplest known cubically convergent process is that due to H. W. 
Richmond [2], viz., 





(2) se iad 26 (en) (xn) 
men" 26" (xn))® — (0) 6" (Xn) 
nati _ 3(¢”(X))? — 26'(X)9""(X) 


It may be shown that for this process lim —— 
. ee ee te! 12(¢"(X))? 





Substituting for ¢’’(x,) in (2) from equation (1), we have 


(3) Xnt1 = Xn — 260" ’ 
2(¢’)*> — os + r¢* + go’ 





where the argument x = x, on the right hand side of the equation (3) is under- 
stood. An illustration of the power of the method is afforded by the computation 
of the first three zeros of Jo(x), which satisfies the differential equation 


xp” + ¢' + xp = 0, 
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ie., p = x,qg = 1,7 =x, s = 0. For this problem ¢’ is easily determined from 
tables since Jo’(x) = — J;(x). Inspection of the tables [3] shows that zeros of 
Jo(x) occur near x = 2.405, x = 5.520, and x = 8.654. Table I gives the results 
obtainable by one application of (3), and contrasts these with those obtained by 
one application of the Newton-Raphson process. 


TABLE I 
Initial Result of Newton-Raphson Result of Richmond's 
approximation process process 
2.405 2.404 825 551 2.404 825 557 697 
5.520 5.520 078 110 5.520 078 110 286 
8.654 8.653 727 909 8.653 727 912 914 


These results were derived solely from tabulated values of Jo(x) and J;(x), and 
required no subsequent interpolation. The quantities in the last column are in 
error in the last figure by one, one and three units, respectively. 

The process is easily extended to the extraction of complex roots. Changing 
to a new independent variable z, where z = x + iy, and letting primes now denote 
differentiation with respect to z, we have, after some manipulation 


' _ DA — BC _ CA + DB 
(4) Yori — Ye = ae Tat — Xn = ae 


where 


A = 2Pre(G re — G im) — 4Pimd' re’ im — GreSre + HimSim + Tre(G re — Pim) 
— 2 imbrim + Ire(breb' re — Dimb' im) — Cim(breb' im + Pim’ re) 
B = 4p, re im + 2Pim(G re — O” im) — PimSre — Grim + 27 rbrbim 
+ Tim(G re — Pim) + Jim(Greh' re — Dim im) + Ire(Greh’ im + Pim re) 
C = 2Dre(dreh' re — Himb' im) + 2Pim(Greh’ im + Pim’ re) 
D = — 2Pim(brb' re — Dimd’ im) — 2Pre(Greb' im + Pim’ re) 


and the real and imaginary parts of p(z,) are written p,, and Pim, and so on. 
An example of the application of the method in the complex case is given by 
the determination of a zero of Yo(z), which satisfies the differential equation 


2p” + o' + 2 = 0; 


ie., we have Pre= Xa, Pim = Yani Gre = 1, Gim =O; Tre = Xap Tim = In; 
Sre = Sim = 0; and Yo'(z) = — Y:(z). Estimation from the tables [4,5] shows 
that a zero of Yo(z) occurs near z = — 2.3762 + 10.6367. (The large number of 
figures in the initial approximation is due to the fact that in these tables Yo(z) 
is tabulated against the variables r and @, where z = re**.) Table II displays the 
results given by the process (4) and contrasts this with the result of one applica- 
tion of the Newton-Raphson process modified for the determination of complex 
roots, viz., 


is or’ re + Dim’ im y -~»y,= ord’ im Pimd’ res 
Ore + o” im in "1 @’ re + o” im 








Xnt+1 — Xn = 
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TABLE II 
Newton-Raphson Richmond’s 
Initial Approximation process process Correct result 
x — 2.3762 — 2.4047 — 2.4031 — 2.4030 
y 0.6367 0.5407 0.5398 0.5399 


It will be noted that this root is not as accurately determined as were the real 
roots of Jo(x) = 0. This is not essentially due to the fact that the computations 
take place in the complex plane, but due to a secondary consequence of this, 
namely that the initial approximation to the root is more inexact. 


P. WYNN 
Scientific Computing Service Ltd. 
23, Bedford Square 
pase W.C. 1, England 
D. R. HartREE, “Notes on iterative processes,” Cambridge Phil. Soc., Proc., v. 45, 1949, 
p. 250-236; see also E. Scuréper, “Uber unendlich viele Algorithmen zir Auflésung der Glei- 
ting Math. Ann., v. 2, 1870, p. 317-363. The possible ambiguity in our definition of con- 
vergence of order k is not si nificant in the present context. 
2. H. W. Ricamonp, “On certain formulae for numerical approximation,’’ London Math. 
ony Jn., v. 19, 1944, p. 31-38. 
. HARVARD University, Tables of the —— Functions of the First Kind of Orders Zero and 
One v. 3, Harvard Univ. Press, Cambridge, 194 
4. CoLuMBIA University, Table of the Bessel Functions Yo(z) and Y,(z) for Complex Argu- 
ments, Columbia Univ. Press, "New York, 195 
5. Cotumpra University, Table of the Bessel Functions Jo(z) and J,(z) for Complex Argu- 
ments, Columbia Univ. Press, 'New York, 1943. 
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33[A].—Epcar Karst, A simple octal (binary) multiplication method for checking 
computer results. 4 leaves, handwritten, deposited in the UMT FILe. 
In this octal multiplication method, sums of one-digit products are accumu- 
lated mentally, with two intermediate lines being written. 
J. L. SELFRIDGE 


University of California 
Los Angeles, California 


34[A-F, HI, K-N, P-V, Z].—Kartw Scniitre, Index Mathematischer Tafelwerke 
und Tabellen (Index of Mathematical Tables), R. Oldenbourg, Miinchen, 1955, 
143, p., 21 cm. Canvas cover. Price, DM 14.50. 


This is a booklet of 143 pages which consists of the following parts: A list of 
approximately 1,200 tables and other useful books classified under 16 general 
headings, an index of authors of these books and tables, and an index of institutes 
referred to in the first part. 

The headings of the 16 classes are as follows: I. Numerical and practical cal- 
culating, II. Logarithms of natural numbers, III. Logarithms of circular func- 
tions, IV. Natural values of circular functions, V. Simple functions derived from 
elementary functions, VI. Primes, prime factors, factors, compound interest and 
rent; continued fractions, theory of numbers, VII. Factorials, gamma functions, 
exponential and hyperbolic functions; elementary transcendental functions, VIII. 
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Elliptic functions and integrals, spherical, Bessel and other higher functions, IX. 
Integral tables and other higher functions; numerical solution of equations and 
differential equations, X. Tables applicable to physics, chemistry, and other 
sciences, XI. Astronomy and astrophysics, XII. Geodesy and geophysics, XIII. 
Nautical and aeronautical determination of position, XIV. Meteorology, XV. 
Astronautics, XVI. Miscellaneous, formulas and tables of measures, weights, 
monetary units. 

There is a preface in both English and German and the Table of Contents is 
also bilingual, as are the headings and subheadings of the first part. In addition, 
there is a list of the abbreviations used in the book. 

The avowed purpose of this volume is essentially the same as that of the earlier 
Index of Mathematical Tables, written by Fletcher, Miller, and Rosenhead, 
(MTAC, 1946-47, v. 2, p. 13-18, 136, 178-181, 219-220, and 277-8). The book 
contains material dated as late as 1955. However, the author states that the 
work was started before World War II and subsequently extended to 1955. Thus, 
it presumably covers 11 years (or almost 11 years) of publication time since the 
publication of Fletcher, Miller, and Rosenhead. 

The 16 large headings are broken down into a total of 130 subdivisions. For 
example, category III. Logarithms of circular functions, has 2 major headings, 
namely, A. Sexagesimal arguments, and B. Centesimals (grades), Sexagesimal and 
centesimal arguments. Under A, the subheadings are degrees, time, degree and 
time, decimals of degree, and decimals of degree and time. Under heading B, the 
subheadings are Centesimal arguments (grades), and Sexagesimal and centesimal 
arguments (grades). 

Within each of the 130 subdivisions (as appropriate), the titles are arranged in 
increasing order of the number of decimal places contained therein and within 
each group where the number of decimal places remains fixed the titles are ar- 
ranged in chronological order. 

Of the approximately 112 pages of indexed tables and books, 42 of them repre- 
sent tables and books in categories 10 through 15, i.e., tables applicable to physics, 
chemistry, biochemistry, biology, paleozoology, technics and engineering, as- 
tronomy and astrophysics, geodesy and geophysics, nautical and aeronautical 
determination of position, meteorology and astronautics. Examination of these 
7 categories leaves the impression that most of the tables indexed hereunder are 
“nonmathematical” in character, i.e., handbooks and special technical tables and 
the like. If the number of titles is approximately proportional to the number of 
pages, this means that only about 750 of the approximately 1200 indexed titles 
are ‘‘mathematical”’ in nature. 

The approximately 450 titles in categories 10-15 are obviously useful informa- 
tion and probably such an index has not heretofore appeared in print. Certainly 
these categories are largely noncompetitive with the earlier index of Fletcher, 
Miller, and Rosenhead. 

With respect to the approximately 750 titles of a “mathematical” nature, a 
sampling shows that probably between 50 and 60 per cent of them are dated 1944 
and earlier and hence are within the purview of the Fletcher, Miller, and Rosen- 
head index. A finer examination of a few categories which can be compared with 
the Fletcher, Miller, Rosenhead index (namely, natural logarithms, natural 
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circular functions (degrees and decimals), hyperbolic functions of real arguments 
and their logarithms, elliptic functions and integrals, Bessel functions of real 
arguments) shows that in the 1944 and earlier category, Schiitte has 42 titles as 
opposed to at least 140 titles in Fletcher, Miller, and Rosenhead. Furthermore, of 
the 42 titles appearing in Schiitte in these categories, 31 also appear in Fletcher, 
Miller, and Rosenhead. If the percentages showing up in this sample persist, this 
would indicate that in the mathematical categories approximately three-fourths 
of the titles which are dated 1944 and earlier are duplicates of those appearing in 
Fletcher, Miller, and Rosenhead. The reviewer had no time to make a complete 
item count and does not know whether the particular categories picked for 
sampling are typical or not. 

The reviewer is somewhat disappointed that many of the desirable features of 
Fletcher, Miller, and Rosenhead index are not included in this book. For example, 
there is no indication whether an item has actually been seen by the author. 
There is no indication of the interval or range of the arguments. There is also 
no indication of the facilities provided for interpolation in tables, or any explana- 
tion of the varying notation, etc. The “standard” or “‘classical’”’ tables in a given 
category are not set forth in boldface type as they are in the F.M.R. index. It 
would seem also that some kind of decimal classification system would be desirable. 

There are some relatively minor things which undoubtedly will be corrected 
in subsequent editions of the book. Some misspellings occur (p. 55, change DeLalla 
to DeLella; p. 83, change Salva to Salvo; p. 138, change J. M. Hammerley to 
J. M. Hammersley). In addition, there are probably a few incorrect titles (for 
example, on p. 55, the book listed as “‘8-figure values of the trigonometrical 
functions at interval 0°.01,”" by E. Buckingham probably is section 1 of the 
Manual of Gear Design written by Mr. Buckingham). A finer breakdown of 
categories would certainly be useful in many instances (for example, in the case 
of hyperbolic functions one might want to know which hyperbolic functions are 
tabulated in a particular place). 

The book will undoubtedly serve a useful purpose since it does have material 
not available elsewhere. However, it is a little disappointing that the same care 
and precision did not go into its preparation as went into the F.M.R. Index, in 
favor of which the reviewer may very well be prejudiced. 

Tuomas H. SouTHARD 
University of California 
Los Angeles, California 


35[E].—NBS Applied Mathematics Series No. 45, Table of Hyperbolic Sines and 
Cosines,x = 2 tox = 10, U.S. Gov. Printing Office, Washington, D. C., 1955, 
v + 81 p., 25.7 cm. Price, 55 cents. 


This table lists sinh x and cosh x for x = 2(.001)10, 9S. 

It is an extension to AMS 36, Tables of Circular Hyperbolic Sines and Cosines 
for Radian Arguments, and its predecessor volume, Mathematical Tables 3, see 
Reviews 89, MTAC, v. 1, 1943, p. 45-7, and 33, MTAC, v. 9, 1955, p. 71-2, and 
note 114, MTAC, v. 4, 1950, p. 123-4. AMS 36 and its predecessors list, among 
other things, sinh x and cosh x for x = 0(.0001)1.9999, 9D. 
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There is adequate introduction. Here the addition formulas used in the 
computation of the table are stated, it is remarked that linear interpolation will 
yield at least six significant digits throughout, and quadratic interpolation formu- 
las accurate to the full nine significant digits are stated. 

In addition, it is remarked in the introduction that $e* is a 9S approximation 
for both sinh x and cosh x for x > 10. 

In the introduction it is claimed that the error in tabulated values should 
be no more than a unit in the last place given. As a matter of routine, Mrs. R. B. 
Horgan and Mr. Edward Sallin recomputed the values for 200 randomly chosen 
arguments taken from [1 ]; they discovered 12 cases in which sinh x disagreed with 
the values in AMS 45 and 6 cases in which cosh x disagreed with the values com- 
puted in AMS 45. Their calculation was based on the values listed in the basic 
table by C. E. Van Orstrand [2]. The discrepancies are all 1 unit in the last place, 
which is within the accuracy claimed in AMS 45. 

The table was computed on the National Bureau of Standards Eastern Auto- 
matic Computer (SEAC) by William F. Cahill. The results are transcribed onto 
punched cards, the values in these cards were then differenced, and the manu- 
script was prepared on the card-controlled typewriter at the National Bureau of 
Standards. The photographic reproduction is excellent, both in format and 
legibility. 

The slight discrepancies noted by Mrs. Horgan and Mr. Sallin are tabulated 
below. 





cS err. 
TABLE 1 
sinh x 
x sinh x from AMS 45 sinh x from SWAC 
2.914 9.18805 728 9.18805 727 
3.388 14.78645 18 14.78645 17 
4.162 32.09210 88 32.09210 87 
4.843 63.42072 31 63.42072 30 
5.404 111.14465 9 111.14465 8 
6.046 211.20880 0 211.20879 9 
6.641 382.92958 3 382.92958 2 
6.882 487.28626 7 487.28626 6 
7.060 582.22215 4 582.22215 3 
8.190 1802.36099 1802.36098 
8.556 2598.92398 2598.92397 
8.572 2640.84121 2640.84120 
TABLE 2 
cosh x 
x cosh x from AMS 45 cosh x from SWAC 
4.617 50.59994 44 50.59994 43 
5.101 82.09605 10 82.09605 09 
5.914 185.09331 8 185.09331 7 
8.510 2482.08165 2482.08164 
8.824 3397.69603 3397.69602 
9,247 5186.69945 5186.69944 
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1. The RAND Corporation, A Million Random Digits with 100,000 Normal Deviates, The Free 
Press, 1955, Review 11, MTAC, v. 10, 1956, p. 39-43. 

2. C. E. VAN OrRsTRAND, “‘Tables of the exponential function and of the circular sine and cosine 
to radian argument,”’ Nat. Acad. Sci. Memoirs, v. 14, 1921, Fifth memoir. 


36[E ].—NBS Applied Mathematics Series No. 46, Table of the Descending Ex- 
ponential, x = 2.5 to x = 10. U. S. Gov. Printing Office, Washington, D. C., 
1955, v + 76 p., 25.7 cm. Price, 50 cents. 


This table lists e~*, x = 2.5(.001)10, 20D. There is a short introduction, which 
describes the computation and printing of the table but otherwise refers the 
reader to an earlier volume [1 ]. 

Calculations were performed on the SEAC computer to 25D, rounded to 20D, 
differenced for checking (no differences are published), and typed by a card-con- 
trolled typewriter (presumably the IBM model 058 at the National Bureau of 
Standards) for photographic reproduction. 

This publication is considered to be auxiliary to [1], which lists (among other 
things) e~*, x = 0(.0001)2.4999, 18D. In the brief introduction the reader is 
referred to this earlier volume for bibliographic notes and other information. 

As a matter of routine, Mrs. R. B. Horgan and Mr. Edward Sallin used the 
SWAC computer to check the values at 200 randomly chosen arguments (taken 
from [2 ]) against products of the values listed by Van Orstrand [3]; there were 
no discrepancies between the listed values and the values computed on SWAC. 

The printing is easily legible and entirely satisfactory. 

. & 3% 


1. NATIONAL BUREAU OF STANDARDS APPLIED MATHEMATICS SERIES 14, Tables of the Exponen- 
tial Function e*, U. S. Government Printing Office, 1951, reviewed in an earlier edition, RMT 524, 
MTAC, v. 3, 1948-9, p. 173. 

2. The RAND Corporation, A Million Random — with 100, 000 Normal Deviates, The Free 
Press, 1955, Review 11, MTAC, v. 10, 1956, p. 39-4 

3. C.E. Van ORsTRAND, ' ‘Tables of the yan function and of the circular sine and cosine 
to radian argument,”’ Nat. Acad. Sci. Memoirs, v. 14, 1921, Fifth memoir. 


37[G].—J. D. Swirt, Tables of Steiner Triple Systems, 4 p., 28 cm., +7 p., 56.cm., 
manuscript. Deposited in UMT FILE. 


These tables list one sytem of each class of isomorphic Steiner triple systems 
of orders 3, 7, 9, 13, and 15. Some algebraic and classification data are included 
in the listings. 

There is only one class for each of the orders 3, 7, and 9; there are two classes 
of order 13; there are 80 classes of order 15. 

The 4 small pages contain introductory and description material and the 
tables through those of order 13. The 7 large pages contain the table of the 80 
classes of order 15. 

Part of the material was computed by the author and a colleague in connection 
with their paper on Steiner triple systems of order 15 [1]. Material was also taken 
from a listing of systems order 13 by F. N. Cole [2]. Cross references are made to 
papers by A. S. White, F. N. Cole, and Louise B. Cummings [3], and R. A. 
Fisher [4]. 

Ci B. T. 


4. phanemeds, HALL HALL, JR., & J. D. Swirt, ‘Determination of Steiner triple systems of order 15,” 
MTAC, v. 9, 1955, p. 146-152. 
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2. F. N. Cog, “The triad systems of 13 letters,” Am. Math. Soc., Trans., 1913, v. 14, p. 1-5. 
3. A. S. Wurtz, F. N. Cote, & Louise D. Cummincs, “Complete classification of etedepeteme 
on . elements,” Nat. Acad. Sci., Memoirs, v. 14, 1925, 2nd Memoir, 89 p. 
. R.A. FISHER, “An examination of the different possible solutions of a problem in incom- 
aa blocks,”” Ann. of Eugenics, v. 10, 1940, p. 52-75. 


38(G, H, I, J, L, M].—Anprew D. Bootn, Numerical Methods, Academic Press, 
New York, Butterworth, London, 1955, viii + 195 p., 83’ & 54”. Price, $6.00. 


Dr. Booth has a reputation [1] in X-ray crystallography, in the design and 
construction of what his publisher [2] calls ‘‘authentic” calculating machines, 
and in their use, e.g., in language translation. Unfortunately this book is not 
up to his usual standards. There are many signs of hasty writing and referencing, 
careless proofreading, and unfamiliarity with current typographical conven- 
tions [3]; indeed Dr. Booth appears to be unsure of the initials of his wife-and 
collaborator (p. 190, 191). In addition there are substantial inaccuracies at various 
levels. 

Among the errors which have caught our eye we mention the following: 

On p. 38, there is the standard classroom mistake of assuming that the mean 
value £, in the error term of a Maclaurin series for f(x), is independent of x, and 
taking it outside the sign of integration. 

On p. 59 and elsewhere we find error terms in approximate quadratures given 
as Ah*f,,, instead of 0Ah"|f"| max with 0 < |@| < 1. 

On p. 92 there is the suggestion that every orthogonal matrix is symmetric. 

We quote, without comment and with the author’s spelling, the first sentence 
of section 7.12 entitled Monte Carlo Methods, and the footnote to which reference 
is made: ‘“The somewhat intriguing title of this section is now applied to methods 
of numerical analysis which make use of the theory of Games [12]. ((12) von 
Neumann, J., and Morganstern, O., The Theory of Games, Princeton (1947)).” 

It is dangerous to recommend methods for the solution of particular problems 
in numerical analysis, such as the inversion of matrices and the determination of 
their eigenvalues, without a knowledge of personnel, equipment available, and of 
the details of the problem. Our experience does not always coincide with that of 
Dr. Booth. 

Nothing less than a complete revision of the book can make it acceptable in its 
purpose to reveal to final honors B.Sc. students at the University of London the 
basic mathematical principles of modern computation. 

The chapter headings are: The Nature and Purpose of Numerical Analysis; 
Tabulation and Differences; Interpolation; Numerical Differentiation and Inte- 
gration; The Summation of Series; The Solution of Ordinary Differential Equa- 
tions; Simultaneous Linear Equations; Partial Differential Equations; Non- 
linear Algebraic Equations; Approximating Functions; Fourier Synthesis and 
Analysis; Integral Equations. 

J. T. 


1. A. D. Bootn, Fourier Technique in X-ray Organic Structure Analysis, Cambridge, 1948. 
A. D. Boor & K. H. V. Boortu, Automatic Digital Calculators, Butterworths Sci. Pub. , London, 
1953. 

W. N. Locke & A. D. Boots, editors, Machine Translation of Language, John Wiley and Sons, 


New York, 1955. 
2. Nature, Butterworths Scientific Publications Advertisement, v. 175, 1955, p. ccccxli. 








106 REVIEWS AND DESCRIPTIONS OF TABLES AND BOOKS 


.. he instance, on one opening, p. 62-63, we find three —_ of indicating the tail of an infinite 
, + etc., «++; 0n p. 53 we find "+ etc., 
wag we find O(x*) signifying a polynomial of dugee ’; aioniitien we find 0(10). 
e reviewer would like to include a disclaimer on the part of SEAC to the discovery of the 
Fenn Mersenne prime. That 2-1 is prime was established by SWAC (see MTAC, 6, 1952, 
p. 205), not by SEAC, as stated on p. 2. 


39[G, H].—Proceedings of the Conference on Linear Programming, May, 1954, 
FERRANTI Ltd., London, England, 1955, 92 mimeographed p., 32.3 cm. A 
limited number of copies are available free from Ferranti Ltd., 21 Portland 
Place, London W. 1, England. 


Ferranti Ltd., which has been constructing high-speed automatic computers, 
issued this report of a conference on linear programming which it sponsored on 
May 4, 1954. Of the eight papers presented, two deal with computation, the re- 
mainder with applications. E. M. L. Beale outlines his ‘‘method of leading vari- 
ables,” previously reported in more detail in [1], and presents worksheets of a 
numerical example. D. G. Prinz describes experiences on the Manchester Com- 
puter with the Simplex Method, using the version of the method in which the 
entire tableau is transformed. He remarks that some difficulties due to round off 
occurred when numbers conceptually zero were not actually zero in the machine. 
These were eliminated by modifying the code so that the conceptual zeros were 
never entered. The paper also offers a brief heuristic justification for the fact that, 
characteristically, the simplex method takes fewer iterations than the conceivable 
upper bound (unknown, by the way), but the argument is not very convincing. 


A. J. HoFFMAN 
National Bureau of Standards 


Washington, D. C. 


1. E. M. L. BEALE, “An alternative method for linear programming,” Camb. Phil. Soc., Proc., 
v. 50, 1954, p. 513-523. 


40[H ].—Marx Lorxin, “A set of test matrices,” 10 p., 27 cm. Typescript de- 
posited in the UMT FIze. 

Table 1 contains inverses of A,, N = 1(1)10. This requires 4 pages. Table 2 
lists the matrices H,, N = 1(1)10. This requires 5 pages. Table 3 lists charac- 
teristic polynomials F,,(A), WN = 1(1)6. This requires about a third of a page. 

The matrix A, is defined in the author’s paper [1], p. 153; H, is defined in 
[1], p. 158; the characteristic polynomials are for the matrices A ,. 


oes. 
1. Marx Lorxin, “A set of test matrices,” MTAC, v. 9, 1955, p. 153-161. 


41[L].—Harotp D. Larsen, Rinehart Mathematical Tables, Formulas and 
Curves, Enlarged Edition, 1953, reprinted 1956. Rinehart & Co., Inc., New 
York, viii + 280 + 2 p., 21 cm. Price $2.50. 


Seven tables, 28-34, have been added to the earlier edition (MTAC, v. 3, 
1949, p. 466) as follows: 


T. 28: Values and Common Logarithms of the Gamma Function, 
T'(n), 5D, m = 1.01(.01)2.50. 


t 
T. 29: The Probability Integral, = f e~"dt, 5D, t = 0(.001)2(.01)3.09. 
0 
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T. 30: Values of Bessel Functions Jo(x) and J:(x), 4D, x = 0(.1)14.9; also 
smallest six positive zeros of Jo(x) and J;(x). 
T. 31: Values of the Complete Elliptic Integrals, 


z/2 
K = f (1 — sin® 6 sin? ¢)~*d¢, 
0 


x/2 
E= f (1 — sin? 6 sin? ¢)d¢, 4D, 
0 
6 = 0°(1°)70°(0.5°)80°(0.2°)89°(0.1°)90°. 


T. 32: Values of the Elliptic Integral of the First Kind, 


F(k, ¢) = f * (1 — sin? @ sin® ¢)-14¢, 4D, 
“0 
8 = 5°(5°)90°, ¢ = 1°(1°)90°. 


T. 33: Values of the Elliptic Integral of the Second Kind, 


: 

E = (k, ¢) = f (1 — sin? @ sin* ¢)*d¢, 4D, 
0 

8 = 5°(5°)90°, ¢ = 1°(1°)90°. 


T. 34: Square Root Divisors (for Calculating Square Roots to 5 Significant 
Figures). 


The typography and composition of the tables are good and pleasing to the 
eye of the user. It is certainly appropriate to have tables of such things as the 
Incomplete Elliptic Integrals of the First and Second Kinds available among a 
collection of elementary tables. 

Several relatively minor things may be worthy of correction in the next edition. 
In T. 29 and T. 32-3, it is not necessary to use the same letter for the upper limit 
and the variable of integration. In T. 32-3, the functions actually tabulated are 
F(6,¢) and E(@,@) respectively instead of those advertised, viz., F(k, ¢) and 
E(k, ¢). For completeness and use in interpolation, it would be nice to have in- 
cluded 6 = 0° and ¢ = 0° among the arguments; however, these entries are easily 
computed and further, their inclusion would have used more space and, in the 
case of 6 = 0°, would probably have necessitated an additional page in each 
table. 

T. 34 appears to be of doubtful utility. With the use of Table 4 and “‘visual 
interpolation,”’ one can usually obtain VN to 4S. Following this with one division, 


1 N 
Newton’s method [ees "2 (xx + *)| gives 7S or 8S. However, the reviewer 
x 
is willing to be convinced by the results of a ‘challenge match.” 


Tuomas H. SouTHARD 


University of California 
Los Angeles, California 
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42[L ].—Gertrupve BiancH & Ipa RuopEs, “Table of characteristic values of 
Mathieu’s equation for large values of the parameter,’’ Washington Academy 
of Sciences, Jn., v. 45, 1955, p. 166-196. 


The tabulation presented in this paper complements the NBS tables (MTAC, 
v. 6, 1952, p. 29-30), and completes the tabulation of characteristic values of 
Mathieu’s equation for all orders r < 15. The notations are the same as in the 
NBS tables, except that the auxiliary parameter ¢ has been used in place of 
s = 1/f, and the quantities 


Be,(#) = be-(s) — (27 + 1)/t, Bo,(#) = bo,(s) — (2r — 1)/t 


have been tabulated in place of the characteristic values. 

The chief aim was a numerical description of the characteristic values for 
r = 0(1)15 and s > 100, or ¢ < .1. The intervals have been so chosen that inter- 
polation by Everett’s formula with second and fourth modified central differences 
will give full accuracy; and the requisite differences are given in the tables. All 
entries are given to 8D, and for ¢ > .01 all entries should be correct, but the 
authors are careful in pointing out that their method of checking made it im- 
possible to detect a random error of two units of the last decimal place. For 
t < .01 (s > 10‘), the accuracy is considerably less, and even the seventh decimal 
place cannot be guaranteed. However, in this range, 6D in B,(#) mean at least 
9S in 5,(#). 

The authors have produced a very useful set of tables. It is to be hoped that 
numerical work on Mathieu functions will continue with emphasis, in the future, 
on functions of higher orders. 


i A. Erp& yi 
California Institute of Technology 


Pasadena, California 


43(Q ].—ApmirALTY HyDROGRAPHIC SUPPLIES ESTABLISHMENT, Tables of Com- 
puted Altitude and Azimuth, v. 1: Latitudes 0°-14° North and South. Hydro- 
graphic Publication H.D. 486/1, Hydrographic Dept., Admiralty, London, 
England, 1952, xvi + 351 p., 29.5 cm. Each volume is priced at 30s. Od. and 
copies may be obtained from Admiralty Chart Agents, the Admiralty Hydro- 
graphic Supplies Establishment, Creechbarrow House, Taunton, and from 
any Admiralty Chart and Chronometer depot. 


The U. S. Hydrographic Office Publication No. 214 has been reproduced in 
Tables of Computed Altitude and Azimuth by the British Hydrographic Depart- 
ment. The original H.O. 214 is in nine volumes, each covering 10° of latitude. The 
regrouping of the new edition results in six volumes, each covering 15°. The 
tabular matter is identical in both, except that certain errors have been corrected 
in the British publication. (These are listed, for Vol. 1, on page xvi of the 
Introduction.) 

The tables are designed to solve directly the spherical triangle of astronomy. 
Latitude, declination, and hour angle are the arguments, yielding altitude and 
azimuth of the object under observation. Computed altitudes less than 5° are 
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omitted, and no refraction is included in the tabulated values. Altitudes are 
printed to the nearest tenth of a minute of arc, azimuths to the nearest tenth 
of a degree. 

Latitudes and hour angles are given at intervals of 1°; declinations at various 
intervals, but not less than 30’. The choice of declination values was made ap- 
propriate to the positions of the Sun, Moon, planets, and the principal navigation 
stars. Thus, the interval is uniformly 30’ up to 30° declination. Above this value 
there are places where very large intervals occur, e.g., between 63° and 69° 
declination. This might cause inconvenience in other applications of the tables 
than that for which they were designed. 

As an aid to interpolation for the actual declination and local hour angle two 
auxiliary quantities Ad and Af are given. These represent, respectively, the change 
of altitude in minutes of arc for one minute change in declination, and the change 
in altitude for a one minute change in hour angle. A multiplication table is printed 
on the end papers to facilitate interpolation. 

A complete description is given of the use of the tables, consistent with naviga- 
tion practice in the British Navy. 


F. H. HOLLANDER 
University of California 
Los Angeles, California 


44[Q ].—ApmrraLty HyproGRAPHIc SupPLIES ESTABLISHMENT, Sight Reduction 
Tables for Air Navigation, v. 1: Selected Stars for Epoch 1955.0, 1952, xi + 324 
p., 29.5 cm.; v. 2: Latitudes 0°-39°, Declinations 0°-29°; v. 3: Latitudes 
40°-89°, Declinations 0°-29°, 1953, x + 338 p., 29.5 cm. A.P. 3270. H.M. 
Nautical Almanac Office, Herstmonceux Castle, Sussex, England. Each vol- 
ume is priced at 30s. Od. and copies may be obtained from Admiralty Chart 
Agents, the Admiralty Hydrographic Supplies Establishment, Creechbarrow 
House, Taunton, and from any Admiralty Chart and Chronometer depot. 


These three volumes represent a joint effort of the U. S. Nautical Almanac 
Office, and H. M. Nautica! Almanac Office. They are issued in the United States 
as Hydrographic Office Publication No. 249. 

Volume 1. Values of altitude, to the nearest minute of arc, and true azimuth, 
to the nearest degree are given for six selected stars for every degree of latitude 
and local hour angle of the first point of Aries. (Above latitude 70° the local hour 
angle is given at intervals of 2°.) The stars were chosen to be those best placed 
for observation at any given time and place. The tables are arranged so as to 
provide for a change in the choice of stars every hour, if necessary (every two 
hours, above 70° latitude). As the star positions are affected by proper motion 
of the stars themselves it is planned to issue a new edition of this volume every 
five years. No refraction has been applied to the tabulated altitudes. 

These tables are of limited value outside their intended area of operation 
which is the planning of observations for aircraft flights, and their subsequent 
reduction. (Various auxiliary tables are given on the end papers. These include 
a table to aid interpolation, one for conversion of angles of arc to time, a refrac- 
tion table, and a table of Coriolis corrections.) 
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The British edition is reproduced from the plates of H.O. 249. The tables were 
originally prepared at the U. S. Nautical Almanac Office of the U. S. Naval 
Observatory. The computations were performed on punched cards, and the 
printers’ copy was prepared on a card-operated typewriter at the Naval Observa- 
tory. Photographic offset reproduction was used. 

Volumes 2 and 3. These are permanent tables for the reduction of aeronautical 
observations of the Sun, Moon, and planets. The arguments are declination: 
0°(1°)29°, and latitude: 0°(1°)89°. The third argument is local hour angle, at 
intervals of 1°(2° above latitude 70°). Volume 2 covers latitudes from 0° to 39°, 
volume 3 from 40° to 89°. Altitude and azimuth are given to the nearest minute 
of arc, and nearest degree, respectively. A column of differences of altitude is 
given for horizontal interpolation from one declination to the next higher one. 

An interesting feature is the extension of the tabulated altitudes to negative 
values for the reduction of solar and lunar observations made at the high levels 
of flight of modern aircraft. A table of the Greenwich Hour Angle and Declination 
of the Sun for the years 1953-2000 is given at the back of each volume. Other 
tables, on the end papers, are similar to those of volume 1. 

The method of preparation, using punched cards, was the same as for the first 
volume. Proofreading on all three volumes was facilitated by the comparison of 
cards punched from proof with those used in the preparation of printers’ copy. 

This helps to explain why the only errata listed are for figures not clearly 
printed, and which may be misread. 


F. H. HOLLANDER 
University of California 


Los Angeles, California 


45(Q ].—SamueE.L Herrick, “Tables for rocket and comet orbits,’’ National 
Bureau of Standards AMS No. 20, U. S. Gov. Printing Office, Washington, 
1953, xxiv + 100 p., 26.7 cm. Price $1.75. 


In Dynamical Astronomy, and more recently in rocketry, the solution is 
frequently required of the equation, E — esin E = M (or the corresponding 
equations for parabolic or hyperbolic motion). This is known as Kepler’s Equa- 
tion, and reputedly more than 300 papers have been written on this subject. Now 
comes another. 

The justification of this one lies in the following considerations. When ‘‘e”’ is 
small, one encounters no difficulty in solving by successive approximations, 
using an ordinary table of sines. As “‘e’’ approaches unity, the number of decimals 
required near the front of the table becomes inordinately large in order to reach a 
comparable accuracy, since one essentially uses (E — sin Z). The author has 
therefore chosen (E — sin £) as the argument of his tables, and the respondents 
are (1 — cos Z£), sin Z, and E (with corresponding functions for the other two 
types of motion). The difficulty with the decimal point cannot be overcome, but 
it is rendered less obtrusive by assigning ‘‘ranges’’ to the argument. 

The Introduction is exceedingly lucid and detailed. The author emphasizes 
the relation of his concept to rectilinear motion, in contrast to the usual approach, 
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which is related to parabolic motion. Illustrations are given to match those of 
Bauschinger’s classical astronomical treatise. 

The tables are clearly printed by photo-offset from copy prepared on a 
punched-card controlled typewriter. Presumably there are no typographical 
errors, if this process is properly prosecuted. The actual calculations were carried 
out under the supervision of Dr. Gertrude Blanch, for the most part on an IBM 
602 Multiplier. The interval of the argument varies throughout the table in such 
a fashion that interpolation is rendered uniform with first and modified second 
differences. Unfortunately, it was not considered feasible to print the one hundred 
entries of the Everett second difference coefficients on one page, or even on one 
opening. 

The application to rockets, as the title imples, is not immediately apparent. 
The only times during which rockets are in rectilinear motion they are subjected 
to such strong forces of thrust or drag that the ordinary equations of two-body 
motion do not apply. The recently proposed artificial satellite will most certainly 
have an “e”’ less than 0.1, and even an intercontinental ballistic missile on any 
kind of a reasonable trajectory has a moderate eccentricity. If some future space 
cadet sets his rocket ship into rectilinear motion with respect to the Sun, he is 
surely doomed to perdition in a very real sense. 


PauL HERGET 
Cincinnati Observatory 
Cincinnati, Ohio 


46(S, T].—Beatrice H. Wors ey, “Solutions of a nonlinear differential equa- 
tion arising in the theory of diffusion flames,’”’ MTAC, v. 9, 1955, p. 112-116. 


The paper describes a numerical method of solving the equation 





PO k —_ 
v sane Tete 1) 


with the following two-point boundary conditious 


v=i-—x, at x«=0 
yv—-0 as x. 


The computations were carried out with the aid of FERUT. Graphs and tables of 
v(x) are available from the author for all the following cases: 


1 — x; = 0.01 with K = 0.001, 0.005, 0.01, 0.05, 0.1, 0.5, 1, and 5. 
1— x, = 0.2 with K = 0.001, 0.005, 0.01, 0.05, 0.1, 0.5, 1, and 5. 
1—x, = 04 with K = 0.005, 0.01, 0.05, 0.1, 0.5, 1, 5, and 10. 
1—x, = 0.6 with K = 0.01, 0.05, 0.1, 0.5, 1, 5, and 10. 
1—x, = 0.8 with K = 0.05, 0.1, 0.5, 1, 5, 10, 50, and 100. 
1—x,=0.9 with K = 0.1, 0.5, 1, 5, 10, 50, and 100. 

A. H. T. 
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47[V].—StTeFan BERGMAN, “Tables for the determination of fundamental solu- 
tions of equations in the theory of compressible fluids,” MTAC, v. 9, p. 8-14, 
1955. 


The tables referred to in the title are tables of the functions F,(A), T(A) for 
complex values of the argument \ where 


3 1-T snes | 
d= $log | =F (+4 h 
T = (i — M} 

—~1\3 
b= (74) 

y+i 


ee ea — 2y)M? — 16] 2 yee 
64 (1 — Mm) <aiin 








F, = (—1)*" 





a; = (37 — 1), a =y7-—3; 


y is the ratio of specific heats presumably taken to satisfy y = 7/5. The argu- 
ments and functional values are given to four decimal places. 


A. H. T. 
NOTES 


Raymond Clare Archibald 
1875-1955 


Professor R. C. Archibald of Brown University, founder and editor of Mathe- 
matical Tables and Other Aids to Computation, died in Sackville, New Brunswick, 
July 26, 1955. He would have been 80 years of age on October 7. 

R. C. A., as readers of MTAC knew him since 1943, had retired as editor in 
1949. A portrait and an account of his life and his writings appeared at that time in 
MTAC, v. 4, 1950, p. 1-2. Since his retirement he had kept up an active interest 
in MTAC contributing articles and reviews of tables. His special interest in 
tables began with his studies of mathematics. In 1908 he commenced the collection 
of tables for the Brown University Library which was to become the most im- 
portant collection of its kind. His connection with the National Research Council 
began in 1939 when he assumed the chairmanship of its MT7AC committee. In a 
prophetic move, he saw the need of a periodical to further the exchange of in- 
formation and ideas about tables and the art of computation. His high standards 
of scientific and bibliographic exactitude have been a challenge to the editors 
who have followed in his steps. 

D. H. LEHMER 


University of California 
Berkeley, California 
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NOTES 


Electronic Computer for the Indian Statistical Institute 


The Hindusthan Standard on November 8, 1955 announced that the Indian 
Statistical Institute will receive an electronic computer from USSR. The Indian 
Statistical Institute is directed by Professor P. C. Mahalanobis. 

Technical details of the computer were not given, but earlier reports by 
members of UNESCO committees indicated that USSR planned to deliver to 
India a copy of the machine built for the Academy of Sciences of the USSR. 
No report of this machine is known to exist in published scientific literature. 

The USSR Academy of Sciences machine was at least partially described at 
the International Conference on Electronic Digital Computers and Information 
Processing held at Darmstadt, Germany, October 25-27, 1955 (announcement, 
MTAC, v. 9, 1955, p. 137). It has also been described in Pravda. 

As it was described in Darmstadt and elsewhere, the machine seems to be a 
floating binary machine with electrostatic memory. The word length is 39 binary 
digits, it uses three address commands, it has approximately 2" words of electro- 
static storage, and a variety of other storage, some novel. 

C.'@..:t. 


Vychislitel’naifa Matematika i Vychislitel’nafa Tekhnika 


We are glad to welcome a foreign competitor. The Institute for Exact Mechan- 
ics and Computational Technology of the Academy of Sciences of USSR began 
recently to issue (irregularly) a journal, Vychislitel’nata Matematika i Vychis- 


litel'nata Tekhnika, that is, Computational Mathematics and Computational Tech- 
nology. Apparently 5000 copies are printed. The first two numbers have reached 
this country. For the time being we propose to publish here a translation of the 
Table of Contents. In accordance with our policy, we do not intend to review the 
articles in numerical analysis but will review articles on tables and on other aids 
to computation, in particular on digital and analog computers. Reviews of the 
articles on numerical analysis are appearing in Mathematical Reviews. 

The translation of the titles in Collection I has been made available by 
Mathematical Reviews ; those of the titles in Collection II have been made by Ida 
Rhodes of National Bureau of Standards. 

It is appropriate to call attention here to the useful sections on numerical 
and graphical methods and on calculating machines and instruments in the 
Russian reviewing journal, Referativnye Zhurnal- Matematika, which appeared first 
in October, 1953. 

i 
Collection I, 1953 


V. A. Ditkin i L. A. Litisternik, Ob odnom prieme prakticheskogo 
garmonicheskogo analiza na sfere [On a method of practical harmonic 
analysis on the sphere ]. 

L. A. Liusternik, Primenenie kubaturnykh formul k chislennomu 
resheniiu. zadachi Koshi dlia nekotorykh uravenenil matematicheskol 
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fiziki [Application of cubature formulas to the numerical solution of 
Cauchy’s problem for certain equations of mathematical physics ]. 

A. I. Ivanova, Nekotorye sluchai kubaturnoi formuly L. A. Lit- 
sternika dlia pravil’nykh mnogougol’nikov [Certain cases of L. A. 
Liisternik’s cubature formula for regular polygons ]. 

A. A. Abramov, O vlifanii oshibok okrugleniia pri reshenii uravneniia 
Laplasa [On the influence of round-off errors in the solution of LaPlace’s 
equation ]. 

L. A. Litisternik, O skhodimosti pri sluchainykh nachal’nykh 
dannykh i nakoplenii oshibok iterafsionnogo profsessa resheniia sistemy 
algebraicheskikh uravnenii [On the convergence and accumulation of 
errors of an iterative process of solution of a system of algebraic equa- 
tions for random initial data ]. 

M. R. Shura-Bura, O reshenii konechno-raznostnogo uravneniia, 
aproksimiruitishchego zadachu Dirikhle dlia uravnenifa Laplasa, na 
élektricheskikh setkakh [On the solution by electric networks of a finite- 
difference equation approximating the Dirichlet problem for the Laplace 
equation ]. 

V. I. Shestakov, Modelirovanie operafsii ischislenifa predlozhenil 
posredstvom prosteishikh chetyrekhpoliisnykh skhem [Modelling the 
operations of the propositional calculus by means of the simplest four- 
pole networks ]. 

A. I. Vzorova, O reshenii sistemy linelnykh algebraicheskikh 
uravnenil sposobom IU. A. Shreidera [On the solution of a system of 
linear algebraic equations by the method of IU. A. Schreider]. 


Collection IT, 1955 


G. S. Khovanskii, Metodika postroenifa nomogramm s orientiro- 
vannym transparantom [A method for the construction of nomograms 
with an oriented transparency ]. 

L. V. Bochek, O chislennom integrirovanil uravnenii v kompleksno! 
oblasti [On the numerical integration of equations in the complex 
domain ]. 

E. N. Dekanosidze, Nekotorye svoistva funkfsii Lommelia ot dvukh 
peremennykh [Certain properties of Lommel functions in two variables. ] 

L. N. Karmazina, Nekotorye svoistva korne! polinomov YAkobi 
[Certain properties of the roots of Jacobi polynomials ]. 

L. N. Karmazina, Ob odnom sposobe vychislenifa gipergeometri- 
chesko!i funkfsii [Concerning one method for the computation of the 
hypergeometric function ]. 

V. K. Saul’ev, K voprosu resheniia zadachi o sobstvennykh znache- 
nifakh metodom konechnykh raznostei [On the question of solving the 
eigenvalue problem by the method of finite differences ]. 

I. M. Stesin, Ofsenka tochnosti vychislenifa sobstvennykh chisel 
apparatom nepreryvnykh drobei [An estimate of the precision of the 
computation of eigenvalues by the means of continued fractions ]. 
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NOTES 


M. A. Tomson, K voprosu ob issledovanii oblastel orbital’noi 
ustoichivosti [On the question of the investigation of regions of orbital 
stability ]. 

M. A. Tomson, Ob odnom sluchae issledovaniia orbital’nof ustoichi- 
vosti reshenil sistemy differenfsial’nykh uravnenif [On one case in the 
investigation of the orbital stability of solutions of systems of differential 
equations ]. 

L. I. Gutenmakher, G. K. Kuz’minok, L. S. Klabukova, Metod 
resheniia lineinykh algebraicheskikh uravnenil na élektronno-lampovom 
integratore [A method of solving linear algebraic equations on a vacuum- 
tube integrator ]. 
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